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Abstract 

In 1864, Waage and Guldberg formulated the "law of mass action." Since that time, chemists, 
chemical engineers, physicists and mathematicians have amassed a great deal of knowledge on the 
topic. In our view, sufficient understanding has been acquired to warrant a formal mathematical 
consolidation. A major goal of this consolidation is to solidify the mathematical foundations 
of mass action chemistry — to provide precise definitions, elucidate what can now be proved, 
and indicate what is only conjectured. In addition, we believe that the law of mass action is of 
intrinsic mathematical interest and should be made available in a form that might transcend its 
application to chemistry alone. We present the law of mass action in the context of a dynamical 
theory of sets of binomials over the complex numbers. 

1 Introduction 

The study of mass action kinetics dates back at least to 1864, when Waage and Guldberg [7] 
formulated the "law of mass action." Since that time, a great deal of knowledge on the topic 
has been amassed in the form of empirical facts, physical theories and mathematical theorems by 
chemists, chemical engineers, physicists and mathematicians. In recent years, Horn and Jackson f9j, 
and Feinberg [5] have made significant mathematical contributions, and these have guided our work. 

It is our view that a critical mass of knowledge has been obtained, sufficient to warrant a formal 
mathematical consolidation. A major goal of this consolidation is to solidify the mathematical 
foundations of this aspect of chemistry — to provide precise definitions, elucidate what can now be 
proved, and indicate what is only conjectured. In addition, we believe that the law of mass action 
is of intrinsic mathematical interest and should be made available in a form that might transcend 
their application to chemistry alone. 

To make the law of mass action available for consideration by researchers in areas other than 
chemistry, we present mass action kinetics in a new form, which we call event-systems. Our formu- 
lation begins with the observation that systems of chemical reactions can be represented by sets 
of binomials. This gives us an opportunity to extend the law of mass action to arbitrary sets of 
binomials. Once this extension is made, there is no reason to restrict ourselves to binomials with 
real coefficients. Hence, we are led to a dynamical theory of sets of binomials over the complex 
numbers. Possible mathematical applications of this theory include: 

1. Binomials are objects of intrinsic mathematical interest [3]. For example, they occur in the 
study of toric varieties, and hence in string theory. With each set of binomials over the 
complex numbers, we associate a corresponding system of differential equations. Ideally, this 
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dynamical viewpoint will help advance the theory of binomials, and enhance our understand- 
ing of their associated algebraic sets. 

2. When we extend the study of the law of mass action to sets of binomials over the complex 
numbers, we can consider reactions that involve complex rates, complex concentrations, and 
move through complex time. Extending to the complex numbers gives us direct access to 
the powerful theorems of complex analysis. Though this clearly transcends conventional 
chemistry, it may have applications in pure mathematics. 

For example, in ongoing work, we seek to exploit an analogy between number theory and 
chemistry, where atoms are to molecules as primes are to numbers. We associate a distinct 
species with each natural number. Then each multiplication rule m x n = ran is encoded 
by a reaction where the species corresponding to the number m reacts with the species 
corresponding to the number n to form the species corresponding to the number mn. With 
an appropriate choice of specific rates of reactions the resulting event-system has the property 
that the sum of equilibrium concentrations of all species at complex temperature s is the value 
of the Riemann zeta function at s. We hope to pursue this approach to study questions related 
to the distribution of the primes. 

3. Systems of linear differential equations are well understood. In contrast, systems of ordi- 
nary non-linear differential equations can be notoriously intractable. Differential equations 
that arise from event-systems lie somewhere in between — more structured than arbitrary 
non-linear differential equations, but more challenging than linear differential equations. As 
such, they appear to be an important new class for consideration in the theory of ordinary 
differential equations. 

In addition to their use in mathematics, event-systems provide a vehicle by which ideas in 
algebraic geometry may be made readily available to the study of mass action kinetics. As such, 
they may help solidify the foundations of this aspect of chemistry. We expand on this in Section [71 

Part of our motivation for this research comes from the emerging field of nanotechnology. To 
quote from [1], "Self-assembly is the ubiquitous process by which objects autonomously assemble 
into complexes. Nature provides many examples: Atoms react to form molecules. Molecules 
react to form crystals and supramolecules. Cells sometimes coalesce to form organisms. Even 
heavenly bodies self-assemble into astronomical systems. It has been suggested that self-assembly 
will ultimately become an important technology, enabling the fabrication of great quantities of small 
objects such as computer circuits. . . Despite its importance, self-assembly is poorly understood." 
Hopefully, the theory of event-systems is a step towards understanding this important process. 

The paper is organized as follows: 

In Section [21 we present the basic mathematical notations and definitions for the study of 
event-systems. 

In Section [3l and all of the sections that follow, we restrict to finite event-systems. Theorem 13.31 
demonstrates that the stoichiometric coefficients give rise to flow-invariant affine subspaces — 
"conservation classes." 

In Sectional and all of the sections that follow, we restrict to "physical event-systems." Though 
we have defined event-systems over the complex numbers, in this paper we focus on consolidating 
results from the mass action kinetics of reversible chemical reactions. Physical event-systems cap- 
ture the idea that the specific rates of chemical reactions are always positive real numbers. The 
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main result of this section is Theorem 14. 5^ which demonstrates that for physical event-systems, if 
initially all concentrations are non-negative, then they stay non-negative for all future real times so 
long as the solution exists. Further, the concentration of every species whose initial concentration 
is positive, stays positive. 

In Section [5l and all the sections that follow, we restrict to "natural event-systems." Natural 
event-systems capture the concept of detailed balance from chemistry. In Theorem l5.lt we give four 
equivalent characterizations of natural event-systems; in particular, we show that natural event- 
systems are precisely those physical event-systems that have no "energy cycles." In Theorem 15. 6^ 
following Horn and Jackson [9], we show that natural event-systems have associated Lyapunov 
functions. This theorem is reminiscent of the second law of thermodynamics. The main result of 
this section is Theorem 15. 15^ which establishes that for natural event-systems, given non-negative 
initial conditions: 

1. Solutions exist for all forward real times. 

2. Solutions are uniformly bounded in forward real time. 

3. All positive equilibria satisfy detailed balance. 

4. Every conservation class containing a positive point also contains exactly one positive equi- 
librium point. 

5. Every positive equilibrium point is asymptotically stable relative to its conservation class. 

For systems of reversible reactions that satisfy detailed balance, must concentrations approach 
equilibrium? We believe this to be the case, but are unable to prove it. In 1972, an incorrect 
proof was offered [Ql Lemma 4C]. This proof was retracted in 1974 [8]. To the best of our knowl- 
edge, this question in mass action kinetics remains unresolved [14:\ p. 10]. We pose it formally in 
Open Problem [H and consider it the fundamental open question in the field. 

In Section [6l we introduce the notion of "atomic event-systems." As the name suggests, this 
is an attempt to capture mathematically the atomic hypothesis that all species are composed of 
atoms. The main theorem of this section is Theorem 16. H which establishes that for natural, atomic 
event-systems, solutions with positive initial conditions asymptotically approach positive equilibria. 
Hence, Open Problem [T] is resolved in the affirmative for this restricted class of event-systems. 

2 Basic Definitions and Notation 

Before formally defining event-systems, we give a very brief, informal introduction to chemical 
reactions. All reactions are assumed to take place at constant temperature in a well-stirred vessel 
of constant volume. 
Consider 

A + 2B C. 

T 

This chemical equation concerns the reacting species A, B and C. In the forward direction, one 
mole of A combines with two moles of B to form one mole of C. The symbol "a" represents a 
real number greater than zero. It denotes, in appropriate units, the rate of the forward reaction 
when the reaction vessel contains one mole of A and one mole of -B. It is called the specific rate of 
the forward reaction. In the reverse direction, one mole of C decomposes to form one mole of A 
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and two moles of B. The symbol "r" represents the specific rate of the reverse reaction. Chemists 
typically determine specific rates empirically. Though irreversible reactions (those with a = or 
T = 0) have been studied, they will not be considered in this paper. 

Inspired by the law of mass action, we introduce a multiplicative notation for chemical reactions, 
as an alternative to the chemical equation notation. In our notation, each chemical reaction is 
represented by a binomial. Consider the following examples. On the left are chemical equations. 
On the right are the corresponding binomials. 

1/3 1 1 



1/2 



3 2 



1/3 1 1 

X^ ^ Xi + X2 7:^3 ~ 7:^1^2 



1/2 



3 2 



2X1 + 3^6 3X1 + 2X2 ^ (^xfxl - txIxI 

T 

Our notation leads us to view every set of binomials over an arbitrary field F as a formal system 
of reversible reactions with specific rates in F \ {0}. For our present purposes, we will restrict our 
attention to binomials over the complex numbers. With this in mind, we now define our notion of 
event-system. 

Notation 2.1. Let Coo = Un^i ^[Xi, X2, • • • , X^^. A monic monomial of Coo is a product of the 
form Y^^i X^^ where the ej are non-negative integers all but finitely many of which arc zero. We 
will write Mqo to denote the set of all monic monomials of Ceo- More generally, if S C {Xi, X2, • • • }, 
we let CfS] be the ring of polynomials with indeterminants in S and we let M5 = M.^ fl CfS] (i.e. 
the monic monomials in CfS"]). 

If n € Z>0! P G C[Xi, X2, • • • , Xn], and a = (ai, a2, • • • , a„) G C"" then, as is usual, we will let 
p{a) denote the value of p on argument a. 

Given two monic monomials M = X^' and AT = H^^i Xf' from Mqo, we will say M 

precedes N (and we will write M ^ N) iS M ^ N and for the least i such that ei fi, ei < fi. 

It follows that 1 is a monic monomial of Cqo and that each element of Coo is a C-lincar 
combination of finitely many monic monomials. Wc will be particularly concerned with the set of 
binomials Boo = {crM + tN | o", t € C \ {0} and M, N arc distinct monic monomials of Cqo}- 

Definition 2.2 (Event-system). An event-system £" is a nonempty subset of Bqo- 

If £ is an event-system, its elements will be called "£"-events" or just "events." Note that if 
aM + tN is an event then M ^ N. 

Our map from chemical equations to events is as follows. A chemical equation 



QiXi ^ bjXj goes to: 



1 -u^t -ru^? if n^r^n^7 

i j i j 

or 2. rHx^ -aHxf^ if Il^^^II^r 



J 
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For example: 

1/2 

1/3 

1/3 

1/2 
-1/2 

-1/3 
-1/2 

1/3 
1/2 

X\ + X2 X^ 
1/3 

3X1 + 2X2 2X1 + 3X6 

T 

Note that our order of monomials is arbitrary. Any linear order would do. The order is necessary 
to achieve a one-to-one map from chemical reactions to events. 

Our definition of event-systems allows for an infinite number of reactions, and an infinite number 
of reacting species. Indeed, polymerization reactions are commonplace in nature and, in principle, 
they are capable of creating arbitrarily long polymers (for example, DNA molecules). 

The next definition introduces the notion of systems of reactions for which the number of 
reacting species is finite. 

Definition 2.3 (Finite-dimensional event-system). An event-system £^ is _^TOte-dimensionaZ iff there 
exists an n G Z>o such that £ C C[Xi, X2, • • • , X„]. 

Definition 2.4 (Dimension of event-systems). Let £^ be a finite-dimensional event-system. Then 
the least n such that £ C C[Xi, X2, • • • , X„] is the dimension of £. 

Definition 2.5 (Physical event, Physical event-system). A binomial e € Bqo is a physical event iff 
there exist a, r G R>o and M, N e Moo such that M ~< N and e = aM - tN. An event-system £ 
is physical iff each e G £^ is physical. 

Chemical reaction systems typically have positive real forward and backward rates. Physical 
event-systems generalize this notion. 

Definition 2.6. Let n G Z>o. Let a = (ai, a2, ■ ■ ■ , an) G C". 

1. q; is a non-negative point iff for i = 1, 2, ... , n, G M>o- 

2. q; is a positive point iff for ? = 1, 2, ... , n, G M>o. 

3. a is a z -point iff there exists an i such that = 0. 

In chemistry, a system is said to have achieved detailed balance when it is at a point where the 
net flux of each reaction is zero. Given the corresponding event-system, points of detailed balance 
corresponds to points where each event evaluates to zero, and vice versa. We call such points 
"strong equilibrium points." 



-X2 - -Xi (because X2 -< Xi) 

^X2 - ^Xi 
1 1 

— X3 — — XiXo 
3 2 

rXfXg - aXfxl 
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Definition 2.7 (Strong equilibrium point). Let £" be a finite-dimensional event-system of dimen- 
sion n. Q G C" is a strong S-equilibrium point iff for all e G f , 6(0;) = 0. 

In the language of algebraic geometry, when £ is a finite-dimensional event-system, its corre- 
sponding algebraic set is precisely the set of its strong f-equilibrium points. 

It is widely believed that all "real" chemical reactions achieve detailed balance. We now intro- 
duce natural event-systems, a restriction of finite-dimensional, physical event-systems to those that 
can achieve detailed balance. 

Definition 2.8 (Natural event-system). A finite-dimensional event-system £ is natural iff it is 
physical and there exists a positive strong (^-equilibrium point. 

Our next goal is to introduce atomic event-systems: finite-dimensional event-systems obeying 
the atomic hypothesis that all species are composed of atoms. Towards this goal, we will define 
a graph for each finite-dimensional event-system. The vertices of this graph are the monomials 
from Moo and the edges are determined by the events. If a weight r is assigned to an edge, then 
r represents the energy released when a reaction corresponding to that edge takes place. For the 
purpose of defining atomic event-systems, the reader may ignore the weights; they are included 
here for use elsewhere in the paper (Definition 15. ip . 

Though graphs corresponding to systems of chemical reactions have been defined elsewhere (e.g. 
[5], [3 p. 10]), it is important to note that these definitions do not coincide with ours. 

Definition 2.9 (Event-graph). Let £ he a finite-dimensional event-system. The event-graph Gs = 
{V, E, w) is a weighted, directed multigraph such that: 

1. F = Moo 

2. For ah Mi, M2 G Moo, for ah r G C, 
(Ml, M2) G ^ and r G u; ((Mi, M2)) iff 

there exist e G f and a, r G C and M,N,T e Moo such that e = aM + tN and M ^ and 
either 

(a) Ml = TM and M2 = TN and r = In (-f ) or 

(b) Ml = TN and M2 = TM and r = - In (-f ) 

Notice that two distinct weights vi and r2 could be assigned to a single edge. For example, let 
£ = {X1X2 — 2Xf,X2 — 5X1}. Consider the edge in Gg from the monomial Xf to the monomial 
XiX2. Weight In 2 is assigned to this edge due to the event X1X2 — 2X1, "^i^^ T = 1. Weight In 5 
is also assigned to this edge due to the event X2 — 5 Ai , with T = Xi. 

Definition 2.10. Let £" be a finite-dimensional event-system. For all M G Moo, the connected 
component of M, denoted Cs{M), is the set of all N G Moo such that there is a path in Gg from 
M to N. 

It follows from the definition of "path" that every monomial belongs to its connected component. 

Definition 2.11 (Atomic event-system). Let £" be a finite-dimensional event-system of dimension 
n. Let S = {Ai,A2,--- ,A„}. Let Ae = {X, e S \ Ce{Xi) = {X,}]. £ is atomic iff for all 
M G Ms, C(M) contains a unique monomial in M^^. 
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If £ is atomic then the members of will be called the atoms of 8. It follows from the definition 
that in atomic event-systems, atoms are not decomposable, non-atoms are uniquely decomposable 
into atoms and events preserve atoms. 

Since the set M|Xi,jf2...,x„} is infinite, it is not possible to decide whether £ is atomic by 
exhaustively checking the connected component of every monomial in Mj^^ X2...,x„}- The following 
is sometimes helpful in deciding whether a finite-dimensional event-system is atomic (proof not 
provided) . 

Let 8 be an event-system of dimension n with no event of the form a + tN. Let = 
{Xi I For ah cr,T G C \ {0} and N G Moo : crXi + tN ^ 8}. Then 8 is atomic iff there exist 
Ml G C'siXi) n Mb£,M2 G Ce{X2) n Mb^, . ..,Mne n such that: 

n n n n 

Forallan^r-^n^'^^^' Il^^'=]l^^'■ W 

1=1 1=1 i=l 1=1 

We have shown (proof not provided) that if 8 and are as above, and there exist Mi G 
Ce{Xi) nMB,,M2 G Ce{X2) nMBe,---,Mn G Ce{Xn)nMBe and there exists a]Xl=iX-^ " 
r nr=i -^i^ ^ ^ such that YYi^=i M"' ^ YYi=i > then 8 is not atomic. Hence, to check whether 
an event-system with no event of the form a + tN is atomic, it suffices to examine an arbitrary 
choice of Mi G CsiXi) nMBg,M2 G Cs{X2) nMB£,...,M„ G Cg{Xn) HMb^, if one exists, and 
check whether ([T]) above holds. 

Example 2.1. Let 8 = {X^ - Xf}. Then Bg = {Xi,X2}. Let Mi = Xi and Ma = X2. 
Trivially, Mi,M2 G M^^, Mi G CsiXi) and M2 G CeiX2). Consider the event X^ - Xf. Since 
M| = X2 ^ Xf = Mf, 8 is not atomic. Note that the event X| — Xl does not preserve atoms. 

Example 2.2. Let 8 = {X} - X2,X| - X:i,X2X^ - Xi}. Then Be = {X^.X^}. Let Mi = 
M2 = X|, M3 = X|, M4 = X4, M5 = X5. Clearly these are all in M^^ . X| - X3 G £: implies 
M3 G ^(Xs). XI-X2 G f: implies M2 G ^(^2). Since (Xi, X2X3, X2XI, X|X|) is a path in G^, 
we have Mi G C£(Xi). For the event Xl - X2, we have M| = X| = M2. For the event X| - X3, 
we have M| = X| = M3. For the event X2X3 - Xi, we have M2M3 = = Mi. Therefore, 8 

is atomic. 

Note that it is possible to have an atomic event-system where is the empty set. For example: 

Example 2.3. Let 8 = {I - Xi}. In this case, 5 = {Xi} and M5 is the set {l,Xi,Xf ,Xf , . . . }. 
It is clear that M5 forms a single connected component C in Gg. Hence, X\ is not in Ag, and 
Ae = 0. 1 is the only monomial in M^^. Since 1 is in C, 8 is atomic. 

3 Finite Event-systems 

The study of infinite event-systems is embryonic and appears to be quite challenging. In the rest 
of this paper only finite event-systems (i.e., where the set 8 is finite) will be considered. It is clear 
that all finite event-systems are finite-dimensional. 

Definition 3.1 (Stoichiometric matrix). Let 8 = {ei, 62, • • • , em} be an event-system of dimension 
n. Let i < n and j < m be positive integers. Let ej = aM + tN, where M ~< N. Then 7j^j is 
the number of times Xi divides minus the number of times Xi divides M. The stoichiometric 
matrix F^ of 8 is the m x n matrix of integers F^- = (7j,j)mxn- 
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Example 3.1. Let d = 0.5X| - bOOXiX^Xj. Let £ = {ei}. Then 71,1 = 1, 71,2 = -2, 71,7 = 
and for all other z, 7i,i = 0, hence r^: = ( 1 -2 1 ). 

Definition 3.2. Let £ = {ei, ■ • • , 6^} be a finite event-system of dimension n. Then: 

1. Ps is the column vector {Pi,P2, . . . , Pn)'^ = r|'(ei, 62, ... , e^)'^. 

2. Let a G C". Then ct is an i--equilibrium point ifi' for i = 1,2, . . . ,n : Pi{a) = 0. 

The Pj's arise from the Law of Mass Action in chemistry. For a system of chemical reactions, the 
Pj's are the right-hand sides of the differential equations that describe the concentration kinetics. 
Definition 13.21 extends the Law of Mass Action to arbitrary event-systems, and hence, arbitrary 
sets of binomials. 

It follows from the definition that for finite event-systems, all strong equilibrium points are 
equilibrium points, but the converse need not be true. 

Example 3.2. Let ei = X2 — and 62 = -^^2 — 2Xi. Let <f = {ei, 62}. Then F^- = f | | j and 



^ = ^ ^- Therefore (2,3) is an £^-equilibrium point. Since ei(2,3) = 1, 



Pi\_f2X2- 3X 
P2 

(2,3) is not a strong iS-equilibrium point. 
Example 3.3. Let ei = 6 - X1X2 and 62 = 2X| - 9Xi. Let £ = {ei, 62}. Then F^ = | 

and P^- = ^ ^ = ^ ^ ^ ^1^^ 1^ ^ ^9-^1 ^ _ rj^-^e point (2,3) is a strong equilibrium point 

because ei(2,3) = and 62(2,3) = 0. Since Pi(2,3) = ei(2,3) +62(2,3) = and P2(2,3) = 
61(2,3) — 262(2,3) = 0, the point (2,3) is also an equilibrium point. 

The event-system in Example 13.21 is not natural, whereas the one in Example 13.31 is. In Theo- 
rem [5]7l it is shown that if £" is a finite, natural event-system then all positive f -equilibrium points 
are strong £'-equilibrium points. 

Definition 3.3 (Event-process). Let £" be a finite event-system of dimension n. Let (Pi,P2, . . . ,Pn)'^ 
Pg. Let C C be a non-empty simply-connected open set. Let / = (/i,/2,-"" , fn) where for 
i = 1, 2, . . . , n, /j : C — > C is defined on $7. Then / is an ^'-process on iff for i = 1,2, ... ,n: 

1. // exists on il. 

2. // = P, o / on n. 

Note that £'-processes evolve through complex time, and hence generalize the idea of the time- 
evolution of concentrations in a system of chemical reactions. 

Definition 13.31 immediatelv implies that if / = (/i, /2, • • • , /n) is an f-process on Q, then for 
i = 1,2, ... ,n, fi is holomorphic on Q. In particular, for each i and all a £ Q, there is a power 
series around a that agrees with fi on a disk of non-zero radius. 

Systems of chemical reactions sometimes obey certain conservation laws. For example, they 
may conserve mass, or the total number of each kind of atom. Event-systems also sometimes obey 
conservation laws. 
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Definition 3.4 (Conservation law, Linear conservation law). Let £ he a finite event-system of 
dimension n. A function g : C" ^ C is a conservation law of £ iff 9 is holomorphic on C", 
g{{0, 0, • • • , 0)) = and Vg ■ Ps is identically zero on C". If g' is a conservation law of £ and g is 
linear (i.e. Vc E C, Vet, f3 E C", g{coL + /3) = cg{cx) + g{P)), then g is a linear conservation law of £■ 

The event-system described in Example 13.21 has a linear conservation law g{Xi,X2) = Xi + X2. 
The next theorem shows that conservation laws of £ are dynamical invariants of ^'-processes. 

Theorem 3.1. For all finite event-systems £, for all conservation laws g of £, for all simply- 
connected open sets C C, for all £-processes f on Q, there exists A; € C such that g o f — k is 
identically zero on $7. 

Proof. Let n be the dimension of £. Let (Pi, P2, . . . , Pn)'^ = Pe- For all t G ^7, by Definition l3.3l for 
i = 1, 2, . . . , n, fi{t) and f-{t) are defined. Further, by Definition l3.4l g is holomorphic on C". Hence, 
gof is holomorphic on Q. Therefore, by the chain rule, {gofy{t) = {Vg\f(^t))-{f[{t), f^{t), . . . , f^it)). 
By Definition [331 for all t e Q, {f{{t),f^{t),...,fUt)) = (Pi(/(t)), P2(/(t)), • • • , Pn(/(t))). From 
these, it follows that {g o /)'(t) = {Vg ■ P£-)(/(t)). But by Definition [331 • Ps is identically 
zero. Hence, for all t € 17, {go /)'(t) = 0. In addition, J7 is a simply-connected open set. Therefore, 
by [21 Theorem 11], there exists A; € C such that g o f — k is identically zero on Q. □ 

The next theorem shows a way to derive linear conservation laws of an event-system from its 
stoichiometric matrix. 

Theorem 3.2. Let £ be a finite event-system of dimension n. For all v G kerF^:, v ■ {Xi, ■ ■ ■ 
is a linear conservation law of £. 

Proof. Let F = F^-, then kerF is orthogonal to the image of F^. By the definition of P = P^-, for 
all w G C", P{w) lies in the image of F-^. Hence, for all v G kerF, for all w G C", v ■ P{w) = 0. 
But V is the gradient of v ■ {Xi, • • • , It now follows from Definition 13.41 that v ■ {Xi, • • • , X„) 
is a linear conservation law of f . □ 

Definition 3.5 (Primitive conservation law). Let £" be a finite event-system of dimension n. For 
all V G kerF^:, the linear conservation law v ■ {Xi,X2, ■ ■ ■ ,Xn) is a primitive conservation law. 

We can show (manuscript under preparation) that in physical event-systems all linear conserva- 
tion laws are primitive and, in natural event-systems, all conservation laws arise from the primitive 
ones. 

Definition 3.6 (Conservation class. Positive conservation class). Let £ he a finite event-system of 
dimension n. A coset of (kerF^-)-*- is a conservation class of £. If a conservation class of £ contains 
a positive point, then the class is a positive conservation class of £. 

Equivalently, a, (3 G C" are in the same conservation class if and only if they agree on all 
primitive conservation laws. Note that if is a conservation class of £ then it is closed in C". The 
following theorem shows that the name "conservation class" is appropriate. 

Theorem 3.3. Let £ be a finite event-system. Let 0, C C be a simply-connected open set containing 
0. Let f be an £-process on Let H be a conservation class of £ containing /(O). Then for all 
t G n, f{t) G H. 
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Proof. Let £, 0, /, H and t be as in the statement of this theorem. For all v G kerF^;, the 
primitive conservation law v ■ {Xi,X2, ■ ■ ■ ,Xn) is a dynamical invariant of /, from Theorem 13.21 
and Theorem 13.11 Hence, 



^'•(/i(0),/2(0),--- 



/n(0)) 



v-{fi{t)j2{t)r-- 



fn{t)) 




•• , fn{0) - fn{t)) = 
fit) G H. 



□ 



4 Finite Physical Event-systems 

In this section, we investigate finite, physical event-systems — a generalization of systems of chem- 
ical reactions. 

It is widely believed that systems of chemical reactions that begin with positive (respectively, 
non-negative) concentrations will have positive (respectively, non-negative) concentrations at all 
future times. This property has been addressed mathematically in numerous papers [H p. 6],[5l 
Remark 3.4], [3, Theorem 3.2], [14, Lemma 2.1]. The notion of "system of chemical reactions" 
varies between papers. Several papers have provided no proof, incomplete proofs or inadequate 
proofs that this property holds for their systems. Sontag [141 Lemma 2.1] provides a lovely proof 
of this property for the systems he considers — zero deficiency reaction networks with one linkage 
class. We shall prove in Theorem 14.51 that the property holds for finite, physical event-systems. 
Finite, physical event-systems have a large intersection with the systems considered by Sontag, but 
each includes a large class of systems that the other does not. We remark that our methods of 
proof differ from Sontag's, but it is possible that Sontag's proof might be adaptable to our setting. 

Lemma 14.41 and Lemma 14.111 are proved here because they apply to finite, physical event- 
systems. However, they are only invoked in subsequent sections. Lemma 14.41 relates £^-processes 
to solutions of ordinary differential equations over the reals. Lemma 14.111 establishes that if an 
^'-process defined on the positive reals starts at a real, non-negative point, then its w-limit set is 
invariant and contains only real, non-negative points. 

The next lemma shows that if two ^'-processes evaluate to the same real point on a real argument 
then they must agree and be real-valued on an open interval containing that argument. The proof 
exploits the fact that £'-processes are analytic, by considering their power series expansions. 

Lemma 4.1. Let £ be a finite, physical event-system of dimension n, let C C be open and 

simply- connected, let f = (/i,/2, • • • , fn) be an £-process on Q and let g = {gi,g2, ■ ■ ■ ,gn) be an 
£ -process on Vl' . // to G H ri' n M and f{to) G M" and /(to) = 9 (to), then there exists an open 
interval /CM such that to ^ I and for all t G /.■ 



1- f{t)=g{t). 



2. For i = l,2 

j G Z>0; Cj 



. . . ,n : if Y^JLq Cj(-z — ^o)"' the Taylor series expansion of fi at to then for all 



3. f{t) G W. 
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Proof. Let k G Z>o. By Definition 13.31 / and g are vectors of functions analytic at to- For 
i = 1, 2, . . . , n, let ^ be the A;*'^ derivative of /, and let /('^^ = (/f \ /f \ . . . , #'^). Define c/f ^ 
and gf'-'^^ similarly. To prove [H it is enough to show that for i = 1,2, ... ,n, fi and gi have the 
same Taylor series around to. Let Vq = {Xi,X2, ■ ■ ■ ,Xn). Let Vk = iSi,c{V k-i)P e (recall that if 
H = {hi{Xi,X2, . . . ,Xm), h2{Xi,X2, ■ ■ ■ ,Xm), ■ ■ ■ , hn{Xi,X2, . . . , Xm)) is a vector of functions in 
m variables then Jac{H) is the n x m matrix (f^), where i = 1,2, ... ,n and j = 1, 2, . . . , m). Let 

{Vk,i, Vk,2, Vk,n) = Vk. We claim that /(^) = o / on 17 and g(^) = Vk o g on i}' and for 
i = 1,2, . . . ,n, Vk^i € R[Xi, X2, . . . , Xrt\. We prove the claim by induction on k. If A; = 0, the proof 
is immediate. If /c > 1, on 

= {Vk-iof)' 
= (Jac(Vfe_i)o/)/' 
= (Jac(Ffe_i)o/)(P^o/) 
= (Jac(Vfe^i)P£)o/ 

= VkOf 

By a similar argument, we conclude that g^^"^ = Vk o g on Q.' . By the inductive hypothesis, 
Vk-i is a vector of polynomials in X2, . . . , It follows that iac{Vk-i) is an n x n 

matrix of polynomials in M[Xi,X2, . . . ,Xri\. Since 8 is physical, Pg is a vector of polynomials in 
M[Xi, X2, . . . , Therefore, Vk = Jac(Vfc_i)P£ is a vector of polynomials in M[Xi, X2, . . . , X„]. 
This establishes the claim. 

We have proved that f^^"^ = Vk o f on Q and g^^^ = Vk o g on 17'. Since, by assumption, 
to G 17 n 17' and /(to) = g{to), it follows that f'^'^^to) = g^''\to). Therefore, for i = 1,2, ... ,n, fi 
and Qi have the same Taylor series around to- For i = 1, 2, . . . , n, let aj be the radius of convergence 
of the Taylor series of fi around to. Let rf = minjgji 2,...,n} (^i- Define Vg similarly. Let Z7 C 17 n 17' 
be some non-empty open disk centered at to with radius r < min(r^, rg). Since 17 and 17' are open 
sets and to G 17 017', such a disk must exist. Letting / = (to — r, to + r) completes the proof of [TJ 

By assumption, /(to) G M", and we have proved that f^^^ = Vk ° f and Vk is a vector of 
polynomials in R[Xi, X2, . . . , X^]. It follows that /(''^(to) G M". Therefore, for i = 1, 2, . . . , n, 
all coefficients in the Taylor series of fi around to are real. It follows that fi is real valued on /, 
completing the proof of El □ 

The next lemma is a kind of uniqueness result. It shows that if two £'-processes evaluate to the 
same real point at then they must agree and be real-valued on every open interval containing 
where both are defined. The proof uses continuity to extend the result of Lemma |4.1[ 

Lemma 4.2. Let £ be a finite, physical event-system of dimension n, let 17,17' C C he open and 
simply- connected, let f = {fi, f2, . . . , fn) be an £ -process on 17 and let g = {gi,g2, ■ ■ ■ , g-n) be an 
£-process on 17'. If G 17 fi 17' and /(O) G and /(O) = g{0), then for all open intervals 
/ C 17 n 17' n M such that G /, for all t G /, /(t) = g{t) and /(t) G M". 

Proof. Assume there exists an open interval / C 17nl7'nM such that G / and S = {t G / | /(t) 7^ 
g{t) or /(t) R'"} ^ 0. Let Bp = B D R>o and let Bn = B D M<o. Note that B = BpU Bn, 
hence, Bp 7^ or Bn ^ 0. Suppose Bp 7^ and let tp = \ni{Bp). By Lemma [4.11 there 



(Inductive hypothesis) 
(Chain-rule of derivation) 
(/ is an ^'-process) 
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exists an e € M>o such that {—e,e) Ci B = 0. Hence, tp > e > 0. By definition of tp, for all 
t G [0,tp), f{t) = g{t) and f{t) € M". Since / and g are analytic at tp, they are continuous at 
tp. Therefore, f{tp) = g{tp) and f{tp) £ M". By Lemma |4.H there exists an e' G M>o such that 
for all t e {tp — e', tp + e'), f{t) = g{t) and f{t) G M", contradicting tp being the infimum of Bp. 
Therefore, Bp = 0. Using a similar agument, we can prove that B]\f = 0. Therefore, B = 0, and 
for all t G /, f{tp) = g{tp) and f{tp) G M". □ 

The next lemma is a convenient technical result that lets us ignore the choice of origin for the 
time variable. 

Lemma 4.3. Let £ he a finite, physical event-system of dimension n, let C C he open and 
simply connected, let f = (/i, /2, . . . , fn) he an £ -process on and let f = {fi, f2, ■ ■ ■ , fn) he an 
£ -process on il.. Let n G and n G and a G M". Let /CM he an open interval. If 

1. f{u) = f{u) = a and 

2. e I and 

3. for all s & I, ti + s G $7 and n + s G 
then for all t G /, f{u + t) = f{u + t). 

Proof. Suppose f{u) = f{u) = a G M". Let r2„ = {z G C [ u+z G il} and VLu = {z G C | u+z G 9}. 
Let h = {hi,h2, . . . ,hn) where for i = l,2,...,n, hi : 0,^ C is such that for all z G ^u, 
hi{z) = fi{u + z) and let h = {hi, h2, . . . , hn) where for i = 1,2, ... ,n, hi : ^ C is such that 
for all z G i^a, hi{z) = fi{u + z). Since u + z is differentiable on and for i = 1, 2, . . . , n, fi 
is differentiable on fi, it follows that for i = 1,2, ... ,n, hi is differentiable on f2„. Further, for 
i = 1, 2, . . . , n, for all z G O.^, h[{z) = f[{u + z) = Pe{fi{u + z)) = Pg{hi{z)), so h is an £'-process 
on Qu- Similarly, h is an f -process on Qa. Note that G i^uf^^a because u G O and ■u G and that 
h{0) = h{0) = a because f{u) = f{u) = ex. By Lemma 14.2^ for all open intervals / C n n M 
such that G /, for all t G /, h{t) = h{t), so f{u + t) = f{u + t). □ 

Because event-systems are defined over the complex numbers, we have access to results from 
complex analysis. However, there is a considerable body of results regarding ordinary differential 
equations over the reals. Definition 14.11 and Lemma [4.41 establish a relationship between ^'-processes 
and solutions to systems of ordinary differential equations over the reals. 

Definition 4.1 (Real event-process). Let be a finite, physical event-system of dimension n. Let 
(Pi, i-*2i • • ■ , Pn)'^ = Ps- Let 7 C R be an interval. Let h = {hi, /12, . . . , hn) where for i = 1, 2, . . . , n, 
/ij : M ^ M is defined on /. Then /i is a real-£'-process on / iff for i = 1,2, . . . ,n: 

1. h'- exists on /. 

2. h'-= Pioh on /. 

Lemma 4.4 (All real-£'-processes are restrictions of f -processes). Let £ be a finite, physical event- 
system of dimension n. Let I CM. be an interval. Let h = {hi, /i2, . . . , hn) he a real-£ -process on 
I. Then there exist an open, simply- connected C C and an £ -process f on such that: 

Lien 
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2. For alltel : f{t) = h{t). 



Proof. Let P = (Pi, P21 ■ ■ ■ , Pn) = Ps- For i = 1,2, ... ,n, Pi is a polynomial and therefore analytic 
on C". By Cauchy's existence theorem for ordinary differential equations with analytic right-hand 
sides [TTj . for all a € /, there exist a non-empty open disk Da Q C centered at a and functions 
fa,i, fa,2, • • • , fa,n analytic on Da such that for i = 1, 2, . . . , n : 

1- fa,iia) = hi{a) 

2. f^^ exists on Da and for all t £ Da : fa^i{t) = Pi{fa,iit), fa,2{t), ■ ■ ■ , fa,n{t)). That is, = 
(/a,i, fa,2, • • • , fa,n) IS an ^-process on Da. 

Claim: For all a G I, there exists 6a G M>o such that for all t S lri{a — 5a,a + 6a) ■ fai^) = h{t). 
To see this, by LemmaHTH for all a G / there exists f3a G IK>o such that for all t G {a— (3a, a+Pa)<^Da, 
fa{t) G Let la = {a — (3a, o- + (3a) H Da. Note that is a real-£'-process on Iq . By the 

theorem of uniqueness of solutions to differential equations with right-hand sides [12], there 
exists 7a G M>o such that for all t G (a — 7a, a -|- 7a) Ci laCi I, fa(t) = h{t). Clearly, we can choose 
5a G M>o such that (a — 5a,a + 5a) ^ (a — 7a, a -|- 7a) n la. This establishes the claim. 

For all a € I, let 5a G M>o be such that for all t £ I (1 {a — 6a, a + 6a) ■ /a(0 = h{t). Let Da be 
an open disk centered at a of radius 5a. 

Claim: For all 01,02 G /, for all t G Da^ n Da2 ■ /ai(^) = /a2(^)- '^o see this, suppose 
Dai ^ ^0.2 7^ ^- I^^t J = Dai ^ 1^ I^- Since Da^ and Da2 are open disks centered on the real 
line, J is a non-empty open real interval. For all t G J, by the claim above, faj^it) = h{t) and 
f a2^^) ~ f^i't)- Hence, fa^it) = fa2(^)- Since J is a non-empty interval, J contains an accumulation 
point. Since fa^ and fa^ are analytic on Da^ n Da2 and Da^ fl L'a2 is simply connected, for all 
t G Dai f~^Da2 '■ /ai(0 = /a2(0- This establishes the claim. 

Let Q = Uaei ^a- Clearly, / C fi. is a union of open discs, and is therefore open. 

For all t G il, there exists a G / such that t G Da. Since Da is a disk, t and a are path-connected 
in J7. Since / is path-connected, and / C f], it follows that 0, is path-connected. 

To see that Q is simply-connected, consider the function R : [0, 1] x Q — > given by (u, z) i— 
Re(2;) -|- ilm{z){l — u). Observe that R is continuous on [0, 1] x Vi, and for all z G 0: -R(0, z) = z, 
i?(l, il) C il, and for all u G [0, 1], for all z G fl M : i?(n, z) G fi. Therefore, i? is a deformation 
retraction. Note that i2(0,il) = 17 and i?(l,0) C M, and 17 is path-connected together imply 
that i?(l,r2) is a real interval. Hence, -R(1,0) is simply-connected. Since R was a deformation 
retraction, 17 is simply-connected. 

Let / : 17 ^ C" be the unique function such that for all a G /, for all t G Da : f{t) = fait)- By 
the claim above and from the definition of 17, / is well-defined. 

Observe that for all t £ I, 



Claim: / is an £'-process on 17. From the definitions of 17 and /, for all t G 17, there exists 
a G / such that t G Da and for all s G Da, f{s) = /a(s)- Since fa is an (^-process on Da, the claim 



Ht) = ftit) 
= fit) 



(Definition of f^) 

(/ C 17 and definition of /). 



follows. 



□ 
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In Theorem l4.5| we prove that if £" is a finite, physical event-system, then f-processes that begin 
at positive (respectively non-negative) points remain positive (respectively non-negative) through 
all forward real time where they are defined. In fact, Theorem 14. 51 establishes more detail about £- 
processes. In particular, if at some time a species' concentration is positive, then it will be positive 
at subsequent times. 

Theorem 4.5. Let £ he a finite, physical event- system of dimension n, let Q Q C be open and 
simply- connected, and let f = (/i, /2, . . . , /«) be an £-process on Vt. If I (1 ^} f) M>o is connected 
and G / and /(O) is a non-negative point then for k = 1,2, . . . ,n either: 

1. For all t € /, fk{t) = 0, or 

2. For alltein M>o, fk{t) G M>o. 

The proof of Theorem 14.51 is highly technical, and relies on a detailed examination of the vector 
of polynomials Pg. This allows us to show (Lemma 14. 8p that if / = (/i, /2, . . . , fn) is an £'-process 
that at real time to is non-negative, then each fi is "right non-negative." That is, the Taylor series 
expansion of fi around to has real coefficients and the first non-zero coefficient, if any, is positive. 
Further, (Lemma 14. lOp if /j(to) = and its Taylor series expansion has a non-zero coefficient, then 
there exists k such that /fc(to) = and the first derivative of fk with respect to time is positive at 



Definition 4.2. Let n G Z>o and let k € {1,2, ... ,n}. A polynomial / G M[Xi, X2, . . . , X„] is 
non-nullifying with respect to k iff there exist m G N, ci,C2,--- ,Cm G ^>o, Mi, M2, ■ ■ ■ , Mm G 
M{Xi,X2,...,x„} and h G X2, . . . , X„] such that / = J2T=i CiM^ + Xkh. 

Observe that for all k, the polynomial is non-nullifying with respect to k. 

Lemma 4.6. Let £ be a finite, physical event-system of dimension n. Let {Pi, P2, ■ ■ ■ , Pn) = Ps- 
Then, for all i G {1, 2, . . . , n}, Pi is non-nullifying with respect to i. 

Proof. Let m = \£\. Let {'yj^i)mxn = ^e- Since £ is physical, there exist ai, (T2, . . . , dm, ti,T2, . . . ,Tm G 
M>o and Mi, M2, . . . , Mm, Ni, N2, . . . , Nm G Moo such that for j = 1,2, ...,m : Mj -< Nj and 
{aiMi - TiNi,a2M2,T2N2, OmM^ - TmNm} = £■ Let i G {1, 2, . . . , n). 

Prom the definition of Ps, Pi = X^jLi 7i,j('7j-Mj- — TjNj). It is sufficient to prove that for 
j = 1,2, ... ,m : 'yj^i{ajMj — rjNj) is non-nullifying with respect to i. Let j G {1,2,..., m}. If 
7j^j = then ^j^i{ajMj — TjNj) = which is non- nullifying with respect to i. If 'jj^i > then, from 
the definition of Tg, Xi \ Nj and 



which is non-nullifying with respect to i since 'yj,iO'j > 0. Similarly, if 'jj^i < then Xi | Mj and 



which is non-nullifying with respect to i since —Jj,iTj > 0. Hence, Pi is non-nullifying with respect 



-fj,i{ajMj - TjNj) = jjAajMj + Xi -JjaTj^ 




-/j,iiajMj - TjNj) 



i 



to i. 



□ 
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Definition 4.3. Let to £ C, let / : C ^ C be analytic at to and let /(t) = X]fet=o ^^(^ ~ ^o)^ 
Taylor series expansion of / around to- Then 0{f,tQ) is the least k such that Ck ^ 0. If for all k, 
Ck = 0, then 0{f,tQ) = oo. 

Definition 4.4 (Right non- negative). Let to G M, let / : C ^ C be analytic at to and let /(t) = 
^^0 '^k{t — to)^ be the Taylor series expansion of / around to- Then / is RNN at to iff both: 

1. For all /c G N, Cfc e M and 

2. Either 0(/,to) = cx) or co{f,to) ^ ^>o- 

Lemma 4.7. Let to G C. Let f,g:C^C be functions analytic at to- Then: 

1. 0{f ■g,to) = 0{f,to) + 0{g,to). 

2. If toe R and f, g are RNN at to then f ■ g is RNN at to- 
The proof is obvious. 

Lemma 4.8. Let 8 he a finite, physical event-system of dimension n, letQ C be open and simply- 
connected and let f = (/i, /2, . . . , fn) be an £-process on $7. For all to € 17 fl M, if /(to) € M"q then 
for i = 1,2, ... ,n : fi is RNN at to. 

Proof Suppose to S n M and /(to) G M|o- Let P = {Pi, P2, . . . , P„) = Pe. Let C = {i \ fi is not 
RNN at to}. 

For the sake of contradiction, suppose C / 0. Let m = minjgc' 0(/j, to). Let A; G C be such 
that O{fk,to) = m. Let fk{t) = Si^o'^«(* ~ ^0)* be the Taylor series expansion of around to. 
Since £ is physical and to G M and /(to) G IR>0) it follows from Lemma 14.1121 that for all i G N, 
aj G M. Further: 

ao = ai = . . . = Om-i = {O{fk,to) = m.) (2) 

am G M<o {fk is not RNN at to.) (3) 

Since /(to) G M>o and G M<o and oo = fkito), it follows that m > 0. 

Consider f^ = Pk ° f ■ By differentiation, the Taylor series expansion of at to is: 

00 

f'At)=Y,ii + l)a,+i{t-to)\ (4) 

i=0 

From Lemma r4.6| Pj. is non- nullifying. Hence, there exist I G N, &i, 62; • • • > ^/ G IR>o, Afi, -^2) • • • j G 
M{Xi,X2,...,x„} and /i G M[Xi, X2, . . . , X„] such that = Y!j=i bjMj +Xk-h. Then for all t G ^i: 

/ 

f'kit) = Pk o fit) = bjM, o fit) + /,(t) • (/i o fit)) (5) 
i=i 

Since h is a polynomial, /i o / is analytic at to. Therefore, /fc ■ (/i o /) is analytic at to. Let 
^^QCj(t — to)* be the Taylor series expansion of /fc • (/i o /) at to- Similarly, for j = 1,2, ... ,Z, 
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bjMj o / is analytic at to- Let YliZo ~ ^o)* be the Taylor series expansion of bjMj o / at to- 
From (Ilj),([5]), equating Taylor series coefficients, for z = 0, l,...,m — 1: 

{i + l)ai+i = Cj + ^ dj^i (6) 

From Lemma l4.7llt 

0{fk ■ {h o /), to) = O{fk,to) + 0{h o /, to) > 0(A, to) = m 

Hence, 

Co = Ci = ... = Cm-1 = 0. (7) 

From ([2]), ([6]), ([7]), for i = 0, 1, . . . , m - 2: 



d,- i = (8) 



Since m > 0, from ([3]), ([6]), ([7]): 



C^i.m-l = "T-Om G K<0 (9) 

Let io = miiij=i,2,...,«{0(6jMj o /,to)}. From Q, it follows that zo < m — 1. 

Case 1: For j = 1, 2, . . . , / : dj^i^^ G M>o. From the definition of io it follows that X^j=i dj^ig G M>o. 
If io < — 1; this contradicts ([8]). If io = m — 1, this contradicts ([9]). 

Case 2: There exists jo G {1,2,...,/} such that dj^^i^ G M<o. From the definition of io, 0{bjQMjQ,to) = 
io < m — 1. Therefore, for each i such that Xi \ Mjg, O{fi,to) < m — 1. From the definitions of C 
and m, this implies that for each i such that Xi \ Mj^, fi is RNN at to. Since bj^ G K>0) it follows 
that bjgMjg o / is a product of RNN functions. Hence, by Lemma 14.7121 bj^AIj^ o / is RNN at to 
and djo,io ^ ^>0^ ^ contradiction. 

Hence, for i = 1, 2, . . . , n, fi is RNN at to. □ 

Lemma 4.9. Let to G M and let f be a function RNN at to- There exists an e € M>o such that 
either for all t G (to, to + e), f{t) G M>o or for all t G (to, to + e), /(t) = 0. 

Proof. Let m = 0{f,t()). If m = oo, / is identically zero and the lemma follows immediately. 
Otherwise, let denote the m*^ derivative of /. Since / is RNN at to and has order m, 

f^"^\to) G M>o. Since / is analytic at to, Z*-™^ is analytic at to, and hence continuous at to- By 
continuity, there exists e G M>o such that for all r G [to, to + e] : f^'^\r) G M>o. From Taylor's 
theorem, for all t G (to, to + e), there exists r G [to, to + such that: 

/(t) = ^^^/'")(r) 

Therefore, /(t) G M>o. □ 
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Note that Lemma 14.81 and Lemma 14.91 together already imply that if £^ is a finite, physical 
event-system, then f-processes that begin at non-negative points remain non-negative through all 
forward real time where they are defined. This result is weaker than Theorem 14.51 

Lemma 4.10. Let £ be a finite, physical event-system of dimension n, let Q Q C be open and 
simply- connected, let f = {fi, /2, . . . , fn) be an £-process on Q.. Let to G If fiio) is non-negative 
and there exists j G {1,2,..., n} such that < O{fj,to) < oo then there exists k E {1, 2, . . . , n} 
such that 0{fk,tQ) = 1. 



Proof Suppose f{to) G M^g. Let C = {i | < O{fi,to) < oo}. Suppose C ^ 0. Let m = 
minjg(7 0(/j, to)- There exists k £ C such that 0{fk,to) = m. 

Let P = {Pi,P2, . . . , Pn) = Pe- From Lemma [4.6l Pk is non-nullifying with respect to k. Hence, 
there exist I en, hiM, ■ ■ ■ M ^ K>o, Mi, M2, ...,Mie M|Xi,x2,...,x„} and h G M[Xi, X2, . . . , X„] 
such that Pk = Yl\=i ^j^j + -^k • h. 

For all t G U: f'^{t) = Pk o f{t) = Y!j=i bjM^ o f{t) + • {h o f{t)). From Lemma I12II1 
0{fk ■ {h o /), to) = 0{fk, to) + 0{h o /, to) > 0{fk, to) = m. It follows that: 



From Lemma (fin2]) and Lemma (jM]), for j = 1,2,... ,/ : bjMj o / is RNN at to. It fol- 
lows that 0{Y^^-^^hjMj o /,to) = m.m.j=i^2,...,iO{hjMj o f,to). From Equation p^ . m — I = 
minj=i^2,. . .,1 0{bjMj o f,to). Hence, there exists jo such that 0{bjgMjQ o /,to) = m — 1. From 
Lemma (|4.7lip . for all i such that Xj | Mj„, 0{fi,tQ) < m — 1. From the definition of m, for all i 
such that Xi \ Mj^, 0{fi,tQ) = 0. It follows that m — 1 = 0(bj^^Mj^^ o /, to) = 0. Hence, m = 1. □ 

We are now ready to prove Theorem 14.51 

Proof of Theorem \4.5\ Suppose / C n M>o is connected and G / and /(O) is a non-negative 
point. If / n M>o = 0, the theorem is immediate. Suppose / fl M>o 7^ 0- 

It is clear that for all k, O(/fc,0) = 00 iff for all t G I, fk{t) = 0. Let C = {i \ O(/i,0) / 00}. 
From Lemma (14. 8p and Lemma (14. 9p . for all k e C, there exists G / fl M>o such that for all 
tG (0,efc) : /fc(t) gM>o. 

Suppose for the sake of contradiction that there exist i G C and t G /nM>o such that /i(t) ^ M>o. 
From Lemma (|4.2p . /j(t) G M. Since fi{ei/2) G M>o and /i(t) G M<o, by continuity there exists 
t' G / n M>o such that fi{t') = 0. 

Let to = inf{t G / fl R>o | There exists i £ C with /j(t) = 0}. It follows that: 

1. to G ]R>o because to > minjgc'{ej}. 
2- /(^o) G M>Q, from the definition of to. 

3. There exists ii G C such that 0(/ji,to) = 1. This follows because there exist io G C 
and T C / n M>o such that to = inf(T) and for all t G T: /jo(t) = 0. By continuity, 
/io(*o) = 0. Hence, 0{fig,to) > 0. Since io G C, O(/j,),0) / 00. By connectedness of /, 
0(/io,to) / 00. Therefore, < O{fi^,to) < 00. Since /(to) G MJ^q, by Lemma (|4.10p . there 
exists ii G {1,2,... ,n} such that 0(/jj,to) = 1. Assume ii ^ C Then O(/j^,0) = 00. By 
connectedness of /, 0(/j^,to) = 00, contradicting that 0(/ii,to) = 1. Hence, ii G C. 




(10) 
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Hence, fi,{to) = 0. Since f{to) G M^q, by Lemma gj]) /^(to) G M>o. 
From the definition of to, for all t G (0, to)) /ii(^) G IR>o- Since to G M>o 




a contradiction. The theorem follows. 



□ 



There is a notion in chemistry that, for systems of chemical reactions, concentrations evolve 
through time to reach equilibrium. In later sections of this paper, we will investigate this notion. 
In the remainder of this section of the paper, we will prepare for that investigation. 

Definition 4.5. Let £" be a finite event-system of dimension n, let 17 C C be open, simply connected 
and such that M>o ^ f2, let / be an f-process on Q, and let q G C". Then q is an oj-limit point 
of f iff for all e G M>o there exists a sequence of non-negative reals {tj}jgz>o such that — > oo as 
i — > oo and for all i G Z>o, ||/(ti) — q\\2 < £■ 

Sometimes, an w-limit is defined by the existence of a single sequence of times such that the 
value approaches the limit. The above definition is easily seen to be equivalent. 

Definition 4.6. Let £^ be a finite event-system of dimension n and let S C C". S is an invariant 
set of £ iff for all q (z S, for all open, simply-connected C C, for all £^-processes / on fi, if G 
and /(O) = q then for ah t G M>o such that [0,t] C n, f{t) G S. 

Lemma 4.11. Let £ be a finite, physical event-system of dimension n, let ^ Q C be open and 
simply connected, and let f be an £-process on Vt. IfM>Q C (7 and /(O) is a non-negative point, 
then the set of all uj-limit points of f is an invariant set of £ and is contained in M>q. 

Proof. Let S be the set of all w-limit points of /. By Lemma l4.5( for all t G M>o, f{t) G M>q, 
hence S C M|g. 

Let q G S", let C C be open, simply-connected, and such that G 12, and let h be an £^-process 
on such that h{0) = q. Suppose u G M>o and [0,u\ C Q. Since £ is finite and physical, Pglm.^ 
can be viewed as a map F : ^ M.-^ of class C^. By Lemma for ah t G [0,u], h{t) G M", so 
^1[o,m] can be viewed as a map X : [0, u] M" such that X' = F{X). By [12^ p. 147], there exists 
a neighborhood U C M" of q and a constant K such that for all a £ U, there exists a unique real-f"- 
process pa defined on [0, u] with pa{0) = ct and \\pa{u) — h{u)\\2 < K\\a — q\\2exp{Ku). Observe 
that necessarily K G M>o. By Lemma 14.41 for all ct G C/ there exists an open, simply-connected 

^ C and an £^-process on such that [0,u] Q and for all t G [0,u], Qait) = Pait)- 
Therefore, ||^ct('u) — h{u)\\2 < K\\a — q\\2ex.];){Ku). 

Let £ G M>o and let 81,62 G M>o be such that K6iex.p{Ku) < e and the open ball centered 
at q of radius 82 is contained in U. Let 6 = min((5i, (52)- Since q is an w-limit point of /, 
there exists a sequence of non- negative reals {ti}jGZ>o such that tj ^ 00 as z — > 00 and for all 
i G Z>o, \\f{ti) — q\\2 < S. Then for all i G Z>o, /(t?) G U, so by Lemma [431 for all t G [0,^], 
f{t, + t) = Qf(^,^){t). Then 



-Fn) - /i(u)||2 = \\Qf^t,){u) - h{u)\\2 



<K\\f{t^)-q\\2exp{Ku) 

< K5 exp{Ku) 

< e 



Thus h{u) is an w-limit point of /, so S is an invariant set of £. 



□ 
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5 Finite Natural Event-systems 



In this section, we focus on finite, natural event-systems — a subclass of finite, physical event- 
systems which has much in common with systems of chemical reactions that obey detailed balance. 

In chemical reactions, the total bond energy of the reactants minus the total bond energy of 
the products is a measure of the heat released. For example, in the reaction, 0X2 — tXi, In (^) is 
taken to be the quantity of heat released. If there are multiple reaction paths that take the same 
reactants to the same products, then the quantity of heat released along each path must be the 
same. 

The finite, physical event-system £ = {2X2 — Xi,X2 — Xi} does not behave like a chemical 
reaction system since, when X2 is converted to Xi by the first reaction. In (2) units of heat are 
released; however, when X2 is converted to Xi by the second reaction, ln(l) = units of heat 
are released. When an event-system admits a pair of paths from the same reactants to the same 
products but with different quantities of heat released, we say that the system has an "energy 
cycle." 

Definition 5.1 (Energy cycle). Let £ he a finite, physical event-system. £ has an energy cycle iff 
G£ has a cycle of non-zero weight. 

Example 5.1. For the physical event-system fi = {2X2 — Xi, X2 — Xi}, the event X2 — X1 induces 
an edge {X2,Xi) in the event graph with weight In (j) = 0. The event 2X2 — Xi induces an edge 
{Xi,X2) with weight — In (y) = — In (2). The weight of the cycle from X2 to Xi and back to X2 
using these two edges, is — In (2) 7^ 0. Hence, £1 has an energy cycle by Definition 15. 1[ 

Example 5.2. For the physical event-system £2 = {X2 — Xi, 2X3X4 — X2X3,X4X5 — X1X5}, 
the cycle (X3X4X5, X2X3X5, X1X3X5, X3X4X5) is induced by the sequence of events 2X3X4 — 
X2X3, X2 — Xi, X4X5 — X1X5 and has corresponding weig ht In f + In i + In i = In (2) / 0. Hence, 
£2 has an energy cycle. 

The following theorem gives multiple characterizations of natural event-systems. 

Theorem 5.1. Let £ be a finite, physical event-system of dimension n. The following are equiv- 
alent: 

1. £ is natural. 

2. £ has a strong equilibrium point that is not a z-point. {i.e. there exists ct G such that for 
all i = 1 to n, 7^ and for all e (z £, e (a) = 0.) 

3. £ has no energy cycles. 

4- If £ = {aiMi — TiXi, CJ2M2 — T2X2, . . . , a^nMrn — TmXm} and for all j = 1 to m, Mj -< Nj 
and crj,Tj > then there exists a £ M" such that T^a = (in (^^^ , . . . , In (^^^ ^ • 

To prove Theorem 15.11 we will use the following lemma. 

Lemma 5.2. Let £ = {aiMi — TiNi,a2M2 — T2N2, . . . ,crmMm — Tf^Nm} be a finite, physical 
event-system of dimension n such that for all j = 1 to m, crj,Tj > and Mj -< Nj. Then for 

alla = (ai,a2,...,a„)^ G M", F^ • a = (in (fj-) , In ),..., In (^) )^ iffie''^--- ,e-") is a 
positive strong £-equilibrium point. 
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Proof. Let £ = {aiMi - TiNi,a2M2 - T2N2, ■■■ , (JmMm - TmNm} and for all j = 1 to m, Mj -< Nj 
and o'j,Tj > 0. Let F = F^-. For all a = {ai, . . . , a„) € M", 



Fq 



In 



,ln 



,ln 



X] '^J'*"* " i^j/'^j) > = 1, 2, . . . , m 



i=l 



, m 



i=l 



^Nj ((e"i , . . . , e"")) /M, ((e°i , . . . , e°")) = c7,-/r,-, Vj = 1, 2, . . . , m 
^ a,Mj ((e"i , . . . , e"")) - r.TV,- ((e^^ , . . . , e"")) = 0, Vj = 1, 2, . . . , m 
44> (e"^ , • • • , e"") is a positive strong £'-equilibrium point. 



Proof of Theorem \5.1\ (4) =^ (1) : Follows from Lemma |5.2[ 

(1) =^ (2) : Follows immediately from definitions. 

(2) ^ (3) : 

Consider an arbitrary cycle C in Gs given by the sequence of k edges 
{{vQ.vi), {vi,V2), . . . , {vk-i,Vk = vo)} with corresponding weights ri,r2, 



(Exponentiation.) 
(Definition of F.) 



□ 



for i = 1, 2, . . . , fc, there exist Tj G Moo and ei ^ £ with Cj 
Mi.Ni G Moo and Mj -< Ni such that either 

1) Vi-i = TiMi and Vi = TiNi and r « = In ^ G tt; ((uj_i, or 

2) Vi-i = TiNi and Vi = TiMi and n = - In ^ G w ((wj-i, fj)) 

Hence, there exists a vector 6 = (61, 62; • • • , ^fc) with 6j = or 1 such that: 



,rfc. By Definition 12. 9^ 
ajMj — TjNj where cr,-,r,- > and 



w 



1=1 1=1 

k k 

{C) = Y,n = Y.{2b,-l)ln 



i=l 



i=l 



(11) 



(12) 



Let a be a strong equilibrium point of £ that is not a z-point. Then, by Definition 12.71 for i = I 

to k, aiMi (q) - TiN, (a) = 

=^ UiMi (q) = TiNi (a) for i = 1 to A; 

^ {aiMi (a))^' = {nNi (a))''' and {nNi = {aiMi {cx))^~'" for i = 1 to A: 

{aiMi {a)t {nN, {cx)f~^^ = {a^Mi («))'"''' {nN, {cx)t for i = 1 to A 
^ n-=i i^iMi {cx)t {nN, {c))'-"^ = nil {cTiMi {cx))'-"^ {TiN, {oc)t 



Wi=l 



bi^.l 



Yli^i ai ^Ti^ [From Equation (1) and since a is not a z-point] 



=^ lli=l ai^-^^Ti^i ~ 

=> X^Li - 1) In = [Taking logarithm] 
^ -u; (C) = [Prom Equation (2)] 
Hence, £ has no energy cycle. 
(3) ^ (4) : 
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Let £ = {(TiMi — TiNi, a2M2 — T2N2, . . . , amMm — "^m-^m} and for all j = 1 to m, Mj -< Nj and 
crj,Tj > 0. Let F = Tg. We shall prove that if the linear equation Ta = (In (ci/ri) , . . . , In {am/Tm))'^ 
has no solution in then £ has an energy cycle. For j = 1 to m, let F-,- be the j*'' row of F. 
If the system of linear equations Fa = (In (cri/ri) , . . . , In {(Tm/'Tm))'^ has no solution in R" then, 
from linear algebra [13^ p. 164, Theorem] and the fact that F is a matrix of integers, it follows that 
there exists Z, there exist (not necessarily distinct) integers ji, ^'2, . . . , G {1, 2, . . . , m}, there exist 
ai, a2, . . . , a; € {+1, —1} such that: 

aiFjj + a2Fj2 H h aiTj^ = (13) 

«i In /-^ii ) + 02 In (crj2 /r,2 ) + •••+ a. In (cJj, /tj, ) 7^ (14) 

Consider the sequence C of Z + 1 vertices in the event-graph defined recursively by 

I I 

vo= n n 

i=l,ai=+l i=l,ai=—l 

and for i = 1 to I, 

7,. — i2_ 

Observe that by (3), 
Hence, 

-0= U u U ^1: U ^" = -^ 

i=l,ai=+l i=l,ai=— 1 j=l,ai=+l i=l,ai=— 1 

Hence, C is a cycle. Further, for z = 1 to 
ttjln^ G ID ((ui.i,^^)) 

From Equation (4), 

w (C) = ai In {(JjjTj^) + 02 In {(JjjTj^) + • • • + In {(JjjTj^) ^ 
Hence, C is an energy cycle. □ 

Horn and Jackson [9j and Feinberg [5j have proved that chemical reaction networks with ap- 
propriate properties admit Lyapunov functions. While finite, natural event-systems are closely 
related to the chemical reaction networks considered by Horn and Jackson and by Feinberg, they 
are not identical. Consequently, we will prove the existence of Lyapunov functions for finite, natural 
event-systems (Theorem 15. 6p . 

The Lyapunov function is analogous in form and properties to "Entropy of the Universe" in ther- 
modynamics. The Lyapunov function composed with an event-process is monotonic with respect 
to time, providing an analogy to the second law of thermodynamics. 
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Definition 5.2. Let £" be a finite, natural event-system of dimension n with positive strong £- 
equilibrium point c = (ci, C2, . . . , c^). Then gs^c '■ ^>o ^ is given by 

n 

gs,c {xi,X2, . . . , x„) = ^ (xj (In {xi) - 1 - In (cj)) + 
i=l 

The function gs^c will turn out to be the desired Lyapunov function. 

Note that if £\ and £2 are two finite natural event-systems of the same dimension and if c is 
a positive strong ^i-equilibrium point as well as a positive strong f2-equilibrium point, then the 
functions gs^^c and g£2,c are identical. 

Lemma 5.3. Let £ = {aiMi ~ tiNi, (T2M2 — T2N2, ■ ■ ■ , amMm — t„iN„i} be a finite, natural event- 
system of dimension n with positive strong £-equilibrium point c, such that for all j = 1 to m, 
aj,Tj > and Mj -< Nj. Then for all x G M"o, 



Vgs,c {x) ■ Ps {x) = f; {a,M, {x) - r,N, (x)) In (^^^^ 
Proof. Let g = gg^c- Let x = {xi,X2, ■ ■ ■ , x„) G K>o- Let P = Pg. 
Vg{x).P{x) = j2{§l:i^)-Pii^)) 

= t.^^{^)(^^jA'^jM,ix)-T,N,ix))^ 
= ^(a,M,(a;)-r,iV,(x))^ln(^(^|j J 

= ^ (a,M, (a;)-r,iV,(a;)) In m(fj j 



i=i \i=i 

TjN, (x) 



J2i<^jMj{x)-TjNj{x))ln 



o-jMj (x) 



The last equality follows from the definition of and the fact that c is a strong-equilibrium 
point. □ 

Lemma 5.4. For all x G M>0; (1 — x) In {x) < with equality iff x = 1. 

Proof If < X < 1 then 1 - a; > and ln(x) < 0. If a; > 1 then 1 - a; < and ln(x) > 0. In either 
case, the product is strictly negative. If a; = 1 then (1 — x) In (x) = □ 

Theorem 5.5. Let £ be a finite, natural event-system of dimension n with positive strong £- 
equilibrium, point c. Then for all x G IR>0' ^9£,c i^) ■ Ps {x) < with equality iff x is a strong 
£-equilibrium point. 
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Proof. Let £ = {(XiMi — tiA''i, 0-2 M2 — T2A''2, . . . , am^m — Tm^m} be a finite, natural event-system 
of dimension n with positive strong ^'-equilibrium point c, such that for all j = 1 to m, (Tj,Tj > 
and Mj -< Nj. Let P = Pg and let g = g^^c- By Lemma 15.31 for all x € M"q, 

Vg (x) ■P{x) = f2 {a,M, (x) - r,N, (x)) In (^^^) 

From Lemma 15.41 and the observation that for j = 1,2, ... ,m, Mj (x) , Nj {x) > when x € IR"q 
and by assumption (Tj,Tj > 0, we have, 

Vg{x) ■Pl.x) <0 

with equality iff for all j = 1,2, ... ,m, (jjA4j{x) = TjNj{x). This occurs iff a; is a strong £- 
equilibrium point. □ 

Recall that a function (7 is a Lyapunov function at a point p for a vector field iff g is smooth, 
positive definite at p and L^g is negative semi-definite at p \1Q\ p. 131]. For a finite natural event- 
system £, Pg induces a vector field on M". We will show that, if c is a positive strong f -equilibrium 
point, then gg^c is a Lyapunov function at c for the vector field induced by Pg. 

Theorem 5.6 (Existence of Lyapunov Function). Let £ he a finite, natural event-system of dimen- 
sion n with positive strong £-equilibrium point c. Then g^^c. is a Lyapunov function for the vector 
field induced by Pg at c. 



Proof. Let g = gs,c- For i = 1,2, . . . ,n: 



dxi V 



which are all in as functions on M"q, hence g \s va. C 

— (c) =ln( — 

dxi \Ci 

establishes that Vg (c) = 0. For z = 1, 2, . . . , n, for A; = 1, 2, . . . , n: 



dxkdxi X 



where 5i^k is the Kronecker delta function. Hence, for all x € M"q, the Hessian of g at x is positive 
definite. Therefore, g is strictly convex over M"q. Further, g (c) = and \7g (c) = and g is strictly 
convex together imply that g is positive definite at c. To establish g as a Lyapunov function, it 
remains to show that the directional derivative Lpg of g in the direction of the vector field induced 
by P = Pg is negative semi-definite at c. This follows from Theorem 15.51 since for all x € R^q, 
Lpg{x) = Vg{x) ■ P{x) <0. □ 

Henceforth, the function gg^c will be called the Lyapunov function of £ at c. The next theorem 
shows that finite, natural event-systems satisfy a form of "detailed balance." 
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Theorem 5.7. If £ is a natural, finite event-system of dimension n then all positive £ -equilibrium 
points are strong £-equilibrium points. 

Proof. Let P = P^. Let c € M"q be a positive strong ^'-equilibrium point. Let ic be a positive 
£'-equilibrium point. That is, P{x) = 0. Hence, Vgg^dx) ■ Ps {x) = 0. By Theorem 15.51 is a 
strong £^-equihbrium point. □ 

The fohowing lemma was proved by Feinberg [5', Proposition B.l]. 

Lemma 5.8. Let n > be an integer. Let U be a linear subspace ofM"^, and let a = (ai, 02, • • • , fln) 
and b be elements 0/ M"o- There is a unique element /j, = (/xi,^2)''' G t^"*" such that 

{aie^^ , 02e'^^ , • • • , CLnG^") — b is an element of U. 

The next theorem follows from one proved by Horn and Jackson [51, Lemma 4B]. Our proof is 
derived from Feinberg's [5, Proposition 5.1]. 

Theorem 5.9. Let £ be a finite, natural event-system of dimension n. Let H be a positive conser- 
vation class of £. Then H contains exactly one positive strong £-equilibrium point. 

Proof. Let P = Pf. Let c* = (cj, C2, . . . , c* ) be a positive strong iS-equilibrium point. Let 
pe Hn RIq. For ah c G RIq, 

(1) c is a strong f -equilibrium point 

^ P(ln(ci), ln(c2), . . . , ln(c„))^ = P(ln(cJ), In(c^), • • • , ln«))^. (LemmaEj]) 

^r(ln(|),ln(|),...,ln(i|))^ = 

44> There exists fi = (/^i, /U2, . . . , G ker P n such that 

(ln(|),ln(|)....,l„(|))" = M. 

<^ There exists /i = {fj-i, fj,2, ■ ■ ■ , fJ-n) ^ ker P n M" such that Cj = c*e^' for i = 1, 2, . . . , n. 

(2) c G F n M" ^ c - p G (ker P)^ n M". (By Definition \3M 

From (1) and (2), c is a positive strong ^'-equilibrium point in H iff there exists /.i G kerP fl M" 
such that c = (c*e'^i , C2e^2^ _ _ _ ^ 0*6^^") and 

(c*e'^i, Cae''^ . . . , c*e^") - p G (kerP)^ n M". Applying Lemma [5^8] with a = c*,b = p and 
U = (kerP)-*- n M", it follows that there exists a unique fi of the desired form. Hence, there 
exists a unique positive strong if-equilibrium point in H given by c = (c^e'^^,C2e'^^, . . . ,c*e'^"). □ 

To prove the main theorem of this section (Theorem l5.15p . we will first establish several technical 
lemmas. 

Lemma 15.101 shows that an event that remains zero at all times along a process can be ignored. 

Lemma 5.10. Let £ be a finite event-system of dimension n, let Q Q C be non-empty, open and 
simply- connected, and let f = (/i, /2, • • • , /n) be an £-process on Q. Then either for all t (z fl, 
f{t) is a strong £-equilibrium point or there exist a finite event-system £ of dimension h < n, an 
£-process / = (/i, /2; • • • > fh) on Q, and a permutation vr on {1, 2, . . . , n} such that: 

1. If £ is physical then £ is physical. 
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2. If S is natural then £ is natural. 



3. If c = {ci, C2, ■ ■ ■ , Cn) is a positive strong £ -equilibrium point, then c = {c^-i (i) , c^-i (2) , • • • , c^-i ( 
is a positive strong £ -equilibrium point. 

4- For all e G there exists t & Q such that e{f{t)) 7^ 0. 

5. If £ is natural, I C $7 n M>o is connected, G / and /(O) is a non-negative point then for all 
t € I n M>o, f{t) is a positive point. 

6. For i = 1,2, . . . ,n, ifTr{i) < h then for all t € Q, fi{t) = f.j^(^i){t). 

7. For i = 1,2, . . . ,n, if irii) > n then for all ti,t2 G ^, fi{ti) = fiit2). 

Proof. Let m = \£\. Let £1 = {e e £ \ there exists t G n,e{f{t)) ^ 0}. U £1 = then for all 
t G r2,e(/(t)) = 0, so f{t) is a strong (£'-equilibrium point and the Lemma holds. Assume f 1 7^ 



and let m = \£i\. For j = l,2,...,m, let (7j,Tj G R>o and Mj = U7=iK'''^^j = DiLi -^i*"'" € 
Moo be such that Mj -< Nj and {aiMi - TiNi,a2M2 - T2N2, . . . ,arhM^ - r^iVrn} = £1 and 

{(TiMi - TiATi, (72M2 - r2iV2, . . . , CFmMm - TmNm} = £■ 

Let C = {i \ there exists j < m such that either aj^i 7^ or bj^i 7^ 0}. Let n = \C\. Let tt be a 
permutation on {1,2,..., n} such that vr(C) = {1,2,..., n}. 

For j = 1,2, . . . ,m, let e^j = aj H^Li ^Z'" '^''-'O' U.i=i -^i''" Let ^ = {e^,i,e^,2, • • • , e,r,m}- 
It follows that £ is a finite event-system of dimension h < n. For i = 1, 2, . . . , n, let /j = f-j^-iiiy 

Let /= (/i, /2, •••,/«)• 

Let (7j,i)mxn = Lf. Let (7j,i)mxn = r^. It follows that for j = 1, 2, . . . , m, for z = 1, 2, . . . , n, 

We claim that / is an ^-process on $7. To see this, for A; = 1, 2, . . . , n, for alH G O : 

f'kit) = f'.-\kp) [Definition of /fc.] 

/ n " \ \ 

[/ is an £^-process on Vt.\ 
[Definition of £^1.] 
[j <m,i^C ^ aj^i = bj^i = 0.] 



m 


-l(fc) 






i=i 




V i=i 


i=l 


m 


-l(fc) 


/ n 


n 


i=i 




\ i=i 


i=l 


m 


-l(fc) 
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i=l 



i=l 



(t) [7r(C7) = {l,2,...,n}.] 



T.^j,n-^{k) ( lliL~Hi)it)P---'^'^ - rj lliL-iii)it))'^---'^^^ ) [By composition. 



i=l 



i=l 



i=l 



i=l 



[From (dH).] 
Definition of /j.] 



j=i 



(t) [Definition of e-i^j- 



Tliis establishes tlie claim. 

With £,h,f and vr as described, we will now establish (1) through (6). 

(1) Follows from the definition of £. 

(2) Follows from 3. 

(3) Suppose £ is natural. Hence, there exists a positive strong ^'-equilibrium point (ci, C2, . . . , c„). 
For j = 1,2, ... ,771 : 



e7r,jK-i(l),C^-l(2), • • • ,C^-i(n)) = 11^^-1 (^) " li^^-MO 

i=l j=l 

iec i£C 

= ej{ci,C2, ... ,Cn) 

= 



[j < 171,1 ^ C ^ aj^i = bj^i = 0.] 



Hence, c is a positive strong f'-equilibrium point. 
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(4) Suppose j < rh. Then for all t € : 

h h 

i=l i=l 
h h 

i=l i=l 

n n 
i=l i=l 

= ej{f{t)) 

Since j < rh, therefore ej € fi and there exists t G J7 such that ej{f{t)) ^ 0. Hence, for all Cttj- € (£", 
there exists i E O such that eTrj{f{t)) / 0. 

(5) Suppose £" is natural, / C J7 n M>o is connected, E I and /(O) is a non-negative point. 
It follows that /(O) is a non-negative point and, from Theorem 14.5^ for all t E /, f{t) is a non- 
negative point. Suppose, for the sake of contradiction, that there exist io < h and Iq G I D IR>o 
such that figito) = 0. From Theorem 14.51 again. fig{0) = and for alH E / : /io(i) = 0. Since I is 
an interval and 0, to £ I, I contains an accumulation point. Hence, since fig is analytic on 0, and 
Q is connected, for all t E O : 

/.oW=0. (16) 

It follows that for all f E : 

rh 

o = fUt) = Y.WnAfit))- (17) 

We claim that for j = 1, 2, . . . , m, for all t £ Q : ^j,ioen,j{f{t)) > 0. 
Case 1: Suppose 'jj^i^ = 0. Then Jj,ioenj{f{t)) = > 0. 
Case 2: Suppose jj^i^ > 0. Then &j.7r-i(io) > 0. Hence, 

n n 

e.,i(/>)) = ^in(/^W)""^""^ -^iH^^W)'"^""^ 

i=l 1=1 
n 

= llikt))"'---'^'^ [Since 6^-, > and fromlM kit) = 0.] 

i=l 

> [f{t) is a non-negative point, by Theorem 14. 5j 
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Hence, 7j,ioe7r,j (fit)) > 0. 

Case 3: Suppose 7j.jp < 0. Then > 0. Hence, 

n n 



,j t \ ~ / / 

i=l 1=1 



-Tj YlifM)''---'^'^ [Since a,- > and from[ISl fi,{t) = 0.] 



1=1 

< [fit) is a non-negative point, by Theorem 14. 5j 

Hence, Jj,ioSTr,j{f{t)) > 0. This completes the proof of the claim. 

From [TTl and the claim, it now follows that for j = 1, 2, . . . , m, for all t G : 

7j,ioe.j(/>)) =0 (18) 

Since io < n, there exists jo < m such that either ajg^ig / or 6jo,io ^ 0. If 7j(,,io ^ then, from 
[13 evrjo(/(0) = 0- If 7jo,jo = th^n, since 7j(,,io = bj^.i^^ - aj^^i^^, it follows that aj^^i^ / and 
bjQ^iQ 7^ 0. Hence, divides Cnja- From [HI it follows that e-K,joif{t)) = 0. Hence, irrespective 
of the value of 7jo,jo' ^'^^ all t € O : enjoifit)) = 0. Since Cj^jg is an element of £, this leads to a 
contradiction with Lemma 15.10141 Hence, for all i < fi, for all t G / fl M>o : fiit) > 0. 

(6) Follows from the definition of /. 

(7) For i = 1,2, ... ,n, if 7r(i) > n then i ^ C. That is, for j = 1,2, ... ,m : jj^i = bj^i — Uj^i = 
— = 0. Hence, for alH G : f'-{t) = X^JLi IjA^jifi^)) = 0. Hence, since fi is analytic on J7, and 
Q is simply-connected, for all ^1,^2 ^ ^ ■ fiih) = /j(^2)- D 

We have described, for finite, natural event-systems, Lyapunov functions on the positive orthant. 
We next extend the definition of these Lyapunov functions to admit values at non- negative points. 

Definition 5.3. Let £" be a finite, natural event-system of dimension n with positive strong £- 
equilibrium point c = {ci,C2, ■ ■ . ,Cn). For all v G M>o, let ^ : M>o ^ M be such that for all 

X G M>o 

fa;(ln(x) - 1 - ln(?j)) -hu, if x > 0; 
9v[x) = < . (19) 

\v, otherwise. 

Then the extended lyapunov function 'gs^ : ]R>q — > M is 

n 

g£^{xi,X2,...,Xn) = '^gr,{xi) (20) 

i=l 

The next lemma lists some properties of extended Lyapunov functions. 

Lemma 5.11. Let £ be a finite, natural event-system of dimension n with positive strong £- 
equilibrium point c = (ci, C2, . . . , c„). Then: 

1- 'gF^ is continuous on M"q. 
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2. For all xi,X2, . . . , 2;„ G M>o, W^i^ii ^2, ■ ■ ■ , Xn) > with equality iff X2, • • • , x„) = c. 

3. For all r £ M>o, the set {x G M>q | 'gs^{x) < r} is bounded. 

4- If ^ C is open, simply connected and such that G il, / = (/i, /2, . . . , fn) is an £-process 
on such that /(O) is a non-negative point, and I C M>o n is an interval such that E / 
then {'g£^ o /) is monotonically non-increasing on I. 

Proof. For i = 1, 2, . . . , n, let be as defined in Equation [T9l 

1. For i = 1,2,... ,n,'g^ is continuous on M>o and lim,j.^Q+ 'g^{x) = Ci = ft^(O), so 'g^ is con- 
tinuous on M>o. Since is the finite sum of continuous functions on M>o, 'gFc is continuous on 

2. Let j € {1, 2, . . . , n}. Let g = 'g^. For all x G M>o, g'{x) = In • If < x < Cj then, 

by substitution, 'g'{x) < 0. Similarly, if x > Cj then > 0. Hence, ^ is monotonically de- 

creasing in (0, Cj) and monotonically increasing in (cj, oo). From continuity of ^ in M>o, it follows 
that 

For all X G M>o,5(x) > g{cj) = with equality iff x = Cj. (21) 
From Equations (j20p and (j2ip . the claim follows. 

3. Observe that \imx—>+oogix) = +oo. It follows that: 

For all r G M>o, the set {x G M>o 1 gix) < r} is bounded. (22) 

xi,X2, . . . ,Xn G M>o are such that 'gsxi{xi, X2, . . . , x„) < r, it follows from Equations (j20l) and 
(j2ip that for i = 1,2,. ..,n : 'g^{xi) < r. The claim now follows from Equation (j22p . 

4. Let C C be open, simply connected, and such that G 0; let / = (/i, /2, • • • , /n) be 
an £^-process on Q such that /(O) is a non-negative point; and let / C M>o fi be an inter- 
val such that G /. By Lemma 15.101 there exists n, £, f, and vr satisfying I5.10llfl5.10l7l Let 
c = (ci, C2, . . . , Cn) = (c^-i(i), c^-i(2), • • • ; Cj^-ij-^)). By Lemma fS. 10121 c is a positive strong equilib- 
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rium point of £. Then for all t G /, 

n 

{W2 o f) it) = Y^gTAfiit)) [Equation m] 

i=l 

i:7r(i)<n i:iT{i)>h 
n 

i=l j:7r(j)>n 
h 

= E^ (•^*(*)) + E ^(Z^^*)) [Definition of c and Lemma EjnEl] 

1=1 i:7r(i)>n 

/)(*)+ E [Equation (1201).] 

«:7r(j)>n 



9e c° f) constant [Lemma lS. 10171 ] 

By Definition [531 for all a; G R^q, 5J7(^e) = 5^ ^(a;). By Lemma I5.1UI51 for alH G / n M>o, 
f{t) G M^o- So for alH G / n M>o, o /) (t) = e ° /) (*)• Then, for alH G / n R>o, 



= V5^,c (/(*)) •/'(*) [Chain rule.] 

= V5£,e (/(*)) • ^£ (/(*)) [Definition ESI] 

<0 [Theorem [531] 



Therefore ° fj non- increasing on / n M>o. 

By Definition 13.31 / is continuous on I; by Theorem 14.51 fil) ^ ^^>0' ^^'^ Lemma 15.11111 
^JT is continuous on M>o; so {wf^° is continuous on /. Therefore {^gj^° is non-increasing 
on /. Thus {'ge^o f) is a constant plus a monotonically non-increasing function on /, so {'gs^o /) 
is monotonically non-increasing on /. □ 

The next lemma makes use of properties of the extended Lyapunov function to show that 
^'-processes starting at non-negative points are uniformly bounded in forward real time. 

Lemma 5.12. Let £ be a finite, natural event-system of dimension n. Let a G M"q. There exists 
k G M>o such that for all^} C C open and simply connected and such that G ri, for all £-processes 
/ = (/i, /2, . . . , fn) on Q such that /(O) = a, for all intervals I Qf] M>o such that G /, for all 
tel, fori = 1,2,..., n: fi{t) G M and < fi{t) < k. 

Proof. Since £ is natural, let c G M"q be a positive strong £'-equilibrium point. Let g = 'gs^. 

Let i = g{a). Let S = {x e M>o | g{x) < £}. By Lemma [5TTT31 S is bounded. Hence, let k be 
such that for all a; G 5 : |a;|oo < k. 
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Let C C be open, simply connected, and such that G ri; let / = (/i, /2, • • • , /n) be an 
£^-process on such that /(O) = a; and let / C IR>o n $7 be an interval such that G /. 
From Theorem 1131 for ah t G /, for i = 1, 2, . . . , n : fi{t) G M and fi{t) > 0. 
Consider the function: 

gof\j:I^R 

From Lemma 15.11141 for all t G /, 5 o /[j is monotonically non- increasing on /. That is, for all 
t G /, 

gifit)) < t (23) 

It follows from Equation [23] and the definition of S that /(/) C S. By the definition of /c, it 
follows that for all t I, for i = 1, 2, . . . , n, fi{t) < k. □ 

The next lemma shows that, because £'-processes starting at non- negative points are uniformly 
bounded in real time, they can be continued forever along forward real time. 

Lemma 5.13 (Existence and uniqueness of £'-process.). Let £ be a finite, natural event-system 
of dimension n. Let a G M>q. There exist a simply- connected open set Q C C, an £ -process 
f = ifii f2, ■ ■ ■ 1 fn) on O and k G M>o such that: 

1. M>o C VL. 

2. /(O) = ex. 

3. For all t G M>o, for i = 1,2, ... ,n : fi{t) andO< fi{t) < k. 

4- For all simply- connected open sets C C, for all £-processes f on Q, for all intervals I C 
Q n M>o, «/ G / and /(O) = a, then for all teL, f{t) = f{t) . 

Proof. Claim: There exists k G M>o such that for all intervals I C M>o with G /, for all real-£^- 
processes h = {hi, /12, . . . , hn) on / with h(0) = a, for all t & I, for i = 1,2, . . . ,n: < hi{t) < k. 

To see this, let / C M>o be an interval such that G /. Let h = {hi, h2, • • • , /in) be a real-f"- 
process on I such that h{0) = a. 

From Lemma|131 there exist an open, simply-connected Q C C and an f -process / = (/i, /2, . . . , /„) 
on 0, such that: 

1. / C f2 

2. For alH G / : f{t) = h{t). 

From Lemma 15.121 there exists k G M>o such that for all t £ I, for i = 1,2, ... ,n: fi{t) G M 
and < fi{t) < k. That is, for all t £ I, for i = 1, 2, . . . , n : < hi{t) < k. This proves the claim. 

Therefore, by HH p. 397, Corollary], there exists k G M>o, there is a real-i5-process h = 
(/ii,/i2, . . . ,hn) on M>o such that h{0) = a and for all t G R>o,for i = 1, 2, . . . , n : < hi{t) < k.. 
By Lemma |4.4| there exist an open, simply-connected C C and an iS-process / on 17 such that 
M>0 C n and for all t G M>o, f{t) = h{t). Therefore, for all t G M>o, for i = 1, 2, . . . , n : fi{t) G M 
and < fi{t) < k. Hence, Parts ()l|2|3p are established. Part([3]) follows from Lemma 14.21 □ 

The next lemma shows that the w-limit points of f-processes that start at non-negative points 
satisfy detailed balance. 



31 



Lemma 5.14. Let £ he a finite, natural event-system of dimension n, letQ (IC be open and simply- 
connected, let f be an £-process on Q, and let q G C". IfM>Q C Q and /(O) is a non-negative 
point and q is an to-limit point of f, then q G M>q and is a strong 8-equilibrium point. 

Proof. Suppose M>o C fi, /(O) is a non-negative point, 5 is the set of w-limit points of /, and 
q € S. By Lemma l4.11l q G K>q. By Lemma 15.131 there exists an open, simply-connected flq Q C 
such that M>o C riq and an £^-process h = {hi, h2, ■ ■ ■ ,hn) on Qq such that h (0) = q. 

Let c be a positive strong £^-equihbrium point. By Lemma l5. 11121 9£,c if (t)) is bounded below 
and, by Lemma 15.11141 is monotonically non-increasing on M>o. Therefore liuit^oo'dFlcif (t)) 
exists. Since is continuous, for all a £ S, 'g£^{a) = limt^oo (/ (i))- By Lemma [4. Ill for 
all t G M>o, h{t) G S. Hence, 'g£^{h{t)) is constant on R>o. 

By Lemma 15.101 either g is a strong £'-equilibrium or there exists a finite event-system £ of 
dimension n < n, an £'-process h = {hi, /12, . . . , h^) on fig, and a permutation vr on {1,2,..., n} 
satisfying 1-7 of Lemma 15.101 

Assume q is not a strong (^-equilibrium point. By Lemma 15.10161 for i = 1,2, . . . ,h, for all t G 
^q, h{t) = /i^-i(i) (t). Let c= {ci,C2,. . . ,Cn) = (c^-i(i),c^-i(2), . . . ,c^-i(n))- By Lemma |5J^ c 
is an ^'-strong equilibrium point. 

For all V G ]R>o, let 'g^ be as defined in Equation [19] in Definition 15. 3[ Then for all t G M>o, 

n h 

&h{t))--gi;^ (h (t)) =Y^g^^{h,{t))-Y,-gE- (hj (t) 

i=i j=i 

n h 

i=i j=i 

n ft 

= Y.9^^)i^--Hi) (*)) - Y.9^^)i^--'U) (*)) 
i=i j=i 

n 

i=fi+l 

But, by Lemma [5. 10171 if tt (i) > n then hi (t) is constant. Hence, 'g^~^ {h^-ij^^-^ (t)) is constant 

for i = n+ 1, n + 2, . . . , n, so 'gs^ {h {t))—'g^ (ji (t)^ is constant. Since {h (t)) and 'gs^ {h (t)) — 

^JT (t)^ are both constant, 'gj~^(ji{t)^ must be constant. By Lemma |5.10I5| for all t G K>o, 

h (t) is a positive point, so by Definitions 15.21 and 15.31 (t)^ = g^ (t)^ . Since g^ - (ji (t)^ 

is constant, ^g^ ^ (ji (t)^ = V^^ ~ (ji {t)^ ■ (ji (t)^ = 0. Then by Theorem 15.51 and continuity 

h{0) must be a strong £'-equilibrium point, so for all e € £, for all t G ilq, e{h{t)) = 0, which 
contradicts Lemma 15.10141 Therefore q is a strong ^'-equilibrium point. □ 

The next theorem consolidates our results concerning natural event-systems. It also establishes 
that positive strong equilibrium points are locally attractive relative to their conservation classes. 
Together with the existence of a Lyapunov function, this implies that positive strong equilibrium 
points are asymptotically stable relative to their conservation classes |10|, Theorem 5.57]. 
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Theorem 5.15. Let 8 he a finite, natural event-system of dimension n. Let H he a positive 
conservation class of £. Then: 

1. For all X G H r\W^Q, there exist k € K>o, an open, simply- connected C C and an £-process 
f = {fi, f2, ■ ■ ■ , fn) on fl such that: 

(a) M>o C n. 

(b) /(O) = X. 

(c) For all t € R>o, f{t) £ Hn W^q. 

(d) For all t G M>o, for i = 1,2, ... ,n, < fi{t) < k. 

(e) For all open, simply- connected J7 C C, for all £-processes f on Q, ifO G $7 and /(O) = x 
then for all intervals I Q Qf] M>o such that G /, for all t (z I : f{t) = f{t). 

2. There exists c £ H such that: 

(a) c is a positive strong £-equilihrium point. 

(h) For all d £ H , if d is a positive strong £-equilihrium point, then d = c. 
(c) There exists U C i7 n M" q such that 
i. U is open in H M" q . 
ii. c G U. 

Hi. For all x , there exist an open, simply- connected C C and an £-process f on 
such that 

A. R>o c n. 

B. /(O) = X. 

C. f (t) ^ c as t ^ oo along the positive real line. (i.e. for all e G M>o, there 
exists to G IR>o such that for all t G M>tQ : \\f{t) — c[|2 < e.) 

Proof. 

1. Follows from Lemma 15.131 and Theorem 13.31 

2. (a) and 2.(b) follow from Theorem 15. 9i 

2.(c) Let c G -fT be a positive strong-£'-equilibrium point as in Theorem 15. 15l2al Let g = gs^c- Let 
T = Hn RIq. For allxeHn M."-, for ah r G M>o, let 

Brix) = {y e H nW \ ll^c - y II2 < r} 
Sr{x) = {y £ H nW \ \\x - y\\2 = r} 
Br{x) = {y € H nW \ \\x - y\\2 < r} 

Since M"q is open in M", it follows that T is open in H D R". Therefore, there exists 6 G M>o 
such that B2s{c) C T. Let 5 G M>o be such that B2s{c) C T. It follows that Bs{c) C T. 

Since g is continuous and Ss{c) is compact, let xq G Ss{c) be such that g{xo) = inix^g^^^f.) oi^)- 
Let U = Bs{c) n {aj G T I g{x) < g{xo)}. It follows that U is open in T. Since xq ^ c, and by 
Lemma [5. 11121 gixo) = 'ge^ixo) > = g{c). Hence, c £ U. 

Let X ^U. From Lemma l5.13( there exist an open, simply-connected 17 C C and an £'-process 
/ on J7 such that IR>o C Q. and /(O) = x. 
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We claim that for all t G M>o, f{t) G Bs{c). Suppose not. Then there exists to £ I^>o such that 
f{to) G Ssic). From the definition of xq, g{xo) < g{f{to)). Since /(O) = x £U, g{f{0)) < g{xo). 
Hence, g{f{0)) < gifito)), contradicting Lemma [5 .11141 

To see that f{t) — > c as t ^ oo along the positive real line, suppose not. Then there exists 
e G M>o such that e < 5 and there exists an increasing sequence of real numbers {tj G M>o}tgz>o 
such that ^ oo as i — 5- oo and for all i, f{ti) G Bs{c) \ B^{c). Since Bs{c) \ Bi,{c) is compact, 
there exists a convergent subsequence. By Definition 14.51 the limit of this subsequence is an w-limit 
point q of / such that q G Bs{c) \ Bi,{c). From Lemma |5.14| g is a strong-f-equilibrium point. 
Since q (z Bs{c), q G T. From Theorem 15.91 q = c. Hence, c ^ -^^(c), a contradiction. □ 

We have established that positive strong equilibrium points are asymptotically stable relative to 
their conservation classes. A stronger result would be that if an (^-process starts at a positive point 
then it asymptotically tends to the positive strong equilibrium point in its conservation class. Such 
a result is related to the widely-held notion that, for systems of chemical reactions, concentrations 
approach equilibrium. We have been unable to prove this result. We will now state it as an open 
problem. This problem has a long history. It appears to have been first suggested in [9| Lemma 4C], 
where it was accompanied by an incorrect proof. The proof was retracted in [8]. 

Open Problem 1. Let 8 he a finite, natural event-system of dimension n. Let H he a positive 
conservation class of £. Then 

1. For all X G H DW^q, there exist k G K>o, an open, simply- connected C C and an £-process 
f = {fi, f2, ■ ■ ■ , fr^ on n such that: 

(a) R>o C n. 
(h) /(O) = X. 

(c) For all t G M>o, f{t) e H n W%. 

(d) For all t G M>o, for i = 1,2, ... ,n, Q < fi{t) < k. 

(e) For all open, simply- connected C C, for all £-processes f on 0,, i/0 G 1^ and /(O) = x 
then for all intervals I Q 0,0 M>o, if £ I then for all t L : f{t) = f{t). 

2. There exists c £ H such that: 

(a) c is a positive strong E-equilihrium point. 

(h) For all d £ H , if d is a positive strong £-equilihrium point, then d = c. 

(c) For all X £ H n M"q, there exist an open, simply- connected 17 C C and an £ -process f 
on such that: 

i. M>o c n. 
a. /(O) = X. 

Hi. f (t) ^ c as t ^ oo along the positive real line. (i.e. for all e G M>o, there exists 
to G M>o such that for all t G ^>to '■ 1 1/(0 ~ c||2 < e.) 

In light of Theorem 15.151 Open Problem [1] is equivalent to the following statement. 

Open Problem 2. Let £ he a finite, natural event-system of dimension n. Let x G IR>o- ^^c'^ there 
exists an open, simply- connected Q Q C, an £-process f on and a positive strong £-equilihrium 
point c such that: 
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1. M>o c n. 



2. /(O) = X. 

3. f (t) ^ c as t ^ oo along the positive real line. (i.e. for all e G IR>0) there exists to ^ IR>o 
such that for all t G M>tg : — c\\2 < e.) 

6 Finite Natural Atomic Event-systems 

In this section, we settle Open [1] in the affirmative for the case of finite, natural, atomic event- 
systems. The atomic hypothesis appears to be a natural assumption to make concerning systems of 
chemical reactions. Therefore, our result may be considered a validation of the notion in chemistry 
that concentrations tend to equilibrium. We will prove the following theorem: 

Theorem 6.1. Let 8 he a finite, natural, atomic event-system of dimension n. Let a G lR>o- 
Then there exists an open, simply- connected il. Q C, an £-process f on il., and a positive strong 
£ -equilibrium point c such that: 

1. M>o C 9., 

2. /(O) = a, and 

3. f{t) — c as t — > CX3 along the positive real line (i.e. for all e G M>o, there exists to G M>o 
such that for all t G ^>to ■ Wfi^) — c|[2 < e). 

It follows from Theorem 15. 151 that the point c depends only on the conservation class of a and 
not on cx itself. That is, two £'-processes starting at positive points in the same conservation class 
asymptotically converge to the same c. 

Implicit in the atomic hypothesis is the idea that atoms are neither created nor destroyed, but 
rather are conserved by chemical reactions. Our proof uses a formal analog of this idea. Recall 
from Definition 12.111 that if £ is atomic then Cs{M) contains a unique monomial from M^^. 

Definition 6.1. Let £" be a finite, natural, atomic event-system of dimension n. The atomic 
decomposition map Dg : M.^Xi,X2,...,Xn} ~^ ^>o function M i— > 62, • • • , &n) such that 

X^'X^' ■■■X^^ £ Cs{M) n ■ 

The next lemma lists some properties of the atomic decomposition map. Note that though the 
event-graph Gs is directed, if M and are monomials and there exists a path in Gs from Af to 

then there also exists a path in Gs from to M. Informally, this is because all events are 
"reversible." 

Lemma 6.2. Let £ he a finite, natural, atomic event-system of dimension n and let M, N G 
M{Xi,X2,...,x„}- Then: 

1. DsiM) = De{N) if and only if Ce{M) = Ce{N). 

2. DeiMN) = Ds{M) + D£{N). 
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Proof. Let D = Ds- 

(1) D{M) = D{N) = (61, 62, ... , bn) if and only ifX^^X^^ • • • ^n" G and X^^X^'' • • • X^" € 
C£(A^). Then Q(M) = Ce{N). 

(2) Let D{M) = (61, 62, • • • , ^n) and D{N) = (ci, C2, . . . , c,i). Then, in G^- there is a path from 
M to X|'iX2^2 . . . x^" G Ma and a path from N to ^2^2 • • • € M^g . It follows that there is a 
path from MN to Xii+'ixf2+''2 . . . jf6n+c„ ^ jjence D{MN) = (61+01,62+02, . . . , 6„+c„) = 
D{M) + D{N). □ 

Definition 6.2. Let £" be a finite, natural, atomic event-system of dimension n. For all i G 
{1,2, . . . for all M G M|Xj^jf2^__jf^}, D£^i{M) is the z*'^ component of De{M). 

Definition 6.3. Let £" be a finite, natural, atomic event-system of dimension n. For all i G 
{1,2,..., n\ the function Ks^i : C" — > C is given by 

n 

(Zl, Z2, ■■■,Zn) ^ > ^ D£^i{Xj)zj. 

i=i 

Lemma 6.3. Xei £ be a finite, natural, atomic event-system of dimension n. Then for all i G 
{1, 2, . . . , n}, the function K^^i is a conservation law of £. 

Proof. Let m = |f |, and for j = 1,2, ... ,m, let o'j,Tj G M>o and A4j, Nj G Mqo with Mj ^ Nj be 
such that <S = {(TiMi — riA^i, . . . , amMm — TmNm}. For z = 1, 2, . . . , n, let aj^i, bj^i G Z>o be such 
that Mj = X^'''X^'-^ ■ ■ ■ Xn'-'' and Nj = X^^^^ X^ ■ ■ ■ X^^'". Let (7J,^)mxn = ^e- 
Then for j = 1,2, ... ,m: 



ajMj - T,Nj G £ 

MjeCeiNj) [Definition ES] 

DeiMj) = Ds{Nj) [Lemma EJ] 

n n 

J2aj,iDe{X,) = Y,hiI^£i^i) [LemmaEj] 

i=l 1=1 

n 

i=l 

n 

J2"fj,i^£i^i) = ^ [Definition [3T] 

i=l 

It follows that for all j G {1, 2, . . . , m}, for all G {1, 2, . . . , n}, 

n 
1=1 

Therefore, for all k G {1,2,..., n}, F^ • {Ds^kiX^), De,k{X2), ■ ■ ■ , De,k{Xn)V = «• Since the 
vector {Dc^ki^i), D£^k{X2), . . ■ ,Ds^k{Xn))'^ is in the kernel of F^-, by Theorem 13.21 Kg^k is a con- 
servation law of □ 



36 



Lemma 6.4. Let £ be a finite, natural event-system of dimension n. Let Af, G Mqo and let 
q G C". If M ^ C£{N) and q is a strong £-equilibrium point and M{q) = 0, then N{q) = 0. 

Proof. Let {vo,vi) be an edge in Gg. Then there exist e £ £ and u, r G M>o and T,U,V G Mqo 
such that e = aU — tV and vq = TU and vi = TV. 

Assume Vo{q) = 0. Then either T{q) = or U{q) = 0. If T{q) = then vi{q) = 0. If C/(q) = 
and q is a strong £'-equihbrium point, then e{q) = aU{q) — TV{q) = 0, so V{q) = 0. Therefore 
vi{q) = 0. The lemma follows by induction. □ 

We are now ready to prove Theorem 16. 1[ 

Proof of Theorem \6.1\ Since a is a positive point, it is in some positive conservation class H. By 
Theorem 15.151 

1. There exists exactly one positive strong if -equilibrium point c £ H. 

2. There exist an open and simply-connected C C and an ^'-process f onO, such that M>o C O 
and /(O) = a. 

3. For ah t G M>o, fit) eHn W^^. 

4. There exists k G M>o such that for i = 1, 2, . . . , n, for ah t G M>o, fi{t) G M and < fi{t) < k. 

Let {tj}j£Zyo be an infinite sequence of non-negative reals such that tj — > oo as j — > oo. Then 
{fitj)}j£Z^o ™ infinite sequence contained in a compact subset of M", so it must have a conver- 
gent subsequence. Let q = {qi, q2, ■ ■ ■ , q-n) G C" be the limit point of a convergent subsequence of 
{/(*i)}iez>o- H and M"q are both closed in C", so g G -ff H M>q. Since £ is natural and q is an 
(j-limit of /, q must be a strong (^-equilibrium point by Lemma 15.141 

Assume, for the sake of contradiction, that q ^ M"q. Let i G {1, 2, . . . , n} be such that qi = 0. 
Let G Cf (Xj) n M^j. Since £ is atomic, a unique such N exists. It follows from the definition 
of event graph that Xi G Ce{N). By Lemma [631 N{q) = Xi{q) = Qi = 0. It follows that N ^ 1. 
Hence, there exists Xa G such that Xa divides N and Xa{q) = 0. 

For all j G {1,2,..., n} such that De,aiXj) / 0, let Mj G Cg^Xj) n M^^. Then Xa divides 
Mj, so Mj{q) = 0. Again by Lemma [Ol Xj{q) = Mj{q) = 0, so qj = 0. It follows that for all 
j G {1, 2, . . . , n} either D£^a{Xj) = or = so 

n 

i=i 

Since K£^a is a conservation law of £ by Lemma 16.31 and q is an w-limit point of /, it follows that 

K£^a{a) = 0. (24) 

For all j, D£^a{Xj) is nonnegative, and o; is a positive point, so for all j G {1, 2, . . . , n}, 
D£^aiXj)aj > 0. But D£^a{Xa) = 1 and > so Ks^aict) > 0, contradicting equation (p^ . 
Therefore g G M"q. Since c is the unique positive strong (^-equilibrium point in H, c = q. 

Let U C f/TlM^Q be the open set stated to exist in Theorem 15. 15l2c1 Since c is an w-limit point 
of /, there exists to ^ I^>o such that /(^o) ^ U. Again by Theorem l5.151 there exist 17 C C and an 
f-process f on Q such that M>o C Q and /(O) = /(to) and /(t) ^ c as t — > oo. By Lemma |4.3| 
for all t G M>o, fit + to) = f{t). Therefore, /(t) ^ c as t ^ oo. □ 
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7 Conclusion 



We have endeavored to place the kinetic theory of chemical reactions on a firm mathematical 
foundation and to make the law of mass action available for purely mathematical consideration. 

With regard to chemistry, we have proven that many of the expectations acquired through 
empirical study are warranted. In particular: 

1. For finite event-systems, the stoichiometric coefficients determine conservation laws that pro- 
cesses must obey (Theorem 13. 3p . In fact, we can show (manuscript in preparation): 

(a) For finite, physical event-systems, the stoichiometric coefficients determine all linear 
conservation laws; 

(b) For finite, natural event-systems, the stoichiometric coefficients determine all conserva- 
tion laws. 

2. For finite, physical event-systems, a process begun with positive (non-negative) concentra- 
tions will retain positive (non-negative) concentrations through forward real time where it 
is defined (Theorem 14. 5p . For finite, natural event-systems, a process begun with positive 
(non-negative) concentrations will retain positive (non-negative) concentrations through all 
forward real time (Theorem 15.150 — that is, it will be defined through all forward real time. 

3. Finite, natural event-systems must obey the "second law of thermodynamics" (Theorem l5.6p . 
In addition, the flow of energy is very restrictive — finite, natural event-systems can contain 
no energy cycles (Theorem 15. ip . 

4. For finite, natural event-systems, every positive conservation class contains exactly one posi- 
tive equilibrium point. This point is a strong equilibrium point and is asymptotically stable 
relative to its conservation class (Theorem 15. 15p . 

Unfortunately, we, like our predecessors, are unable to settle the problem of whether a process 
begun with positive concentrations must approach equilibrium. We consider this the fundamental 
open problem in the field (Open Problem [1]). For finite, natural event-systems that obey a math- 
ematical analogue of the atomic hypothesis, we settle Open Problem [T] in the affirmative (The- 
orem [6Tj. In particular, we show that for finite, natural, atomic event-systems, every positive 
conservation class contains exactly one non-negative equilibrium point. This point is a positive 
strong equilibrium point and is globally stable relative to the intersection of its conservation class 
with the positive orthant. 

In terms of expanding the mathematical aspects of our theory, there are several potentially 
fruitful avenues including: 

1. Complex-analytic aspects of event-systems. While we exploit some of the complex- 
analytic properties of processes in this paper, we believe that a deeper investigation along 
these lines is warranted. For example, if we do not restrict the domain of a process to be 
simply-connected, then each component of a process becomes a complete analytic function in 
the sense of Weierstrass. 

2. Infinite event-systems. Issues of convergence arise when considering infinite event-systems. 
To obtain a satisfactory theory, some constraints may be necessary. For example, a bound on 
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the maximum degree of events may be worth considering. It may also be possible to generalize 
the notion of an atomic event-system to the infinite-dimensional case in such a way that each 
atom has an associated conservation law. One might then restrict initial concentrations to 
those for which each conservation law has a finite value. Additional constraints are likely to 
be needed as well. 

3. Algebraic-geometric aspects of event-systems. Every finite event-system that generates 
a prime ideal has a corresponding affine toric variety (as defined in [U p. 15]). The closed 
points of this variety are the strong equilibria of the event-system. Further, every affine toric 
variety is isomorphic to an affine toric variety whose ideal is generated by a finite event system. 
One could generalize event-systems to allow irreversible reactions. In that case, it appears 
that the prime ideals generated by such event-systems are exactly the ideals corresponding 
to affine toric varieties. 

We can show (proof not provided) that finite, natural, atomic event-systems generate prime 
ideals. We are working towards settling Open Problem [1] in the affirmative for every finite, 
natural event-system that generates a prime ideal. 
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